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METHOD AND APPARATUS FOR IMPROVING RESOLUTION 
IN SPECTROMETERS PROCESSING OUTPUT STEPS FROM 
NON-IDEAL SIGNAL SOURCES 

5 CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application is a continuation-in-part of U.S. Patent Application No. 
10/429,170, filed May 1, 2003 (William K. Warburton and Michael Momayezi), which is 
continuation of U.S. Patent Application No. 09/643,705, filed August 21, 2000 (WiUiam 
10 K. Warburton and Michael Momayezi), now U.S. Patent No. 6,587,814, which claims 
priority from Provisional U.S. Patent Application No. 60/151,020, filed August 27, 1999 
(William K. Warburton and Michael Momayezi), all of whose disclosures are 
incorporated by reference in their entirety for all purposes. 

1 5 STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER 

FEDERALLY SPONSORED RESEARCH AND DEVELOPMENT 

[0002] The U.S. Government has rights in this invention pursuant to Contract No. DE- 
FG03-97ER82510 av/arded by the Department of Energy. 

20 BACKGROUND OF THE INVENTION 

[0003] The present invention relates generally to signal processing systems, and more 
particularly to processing the step-lilce output signals generated by devices with multiple 
poles and zeros (including non-ideal, nominally single-pole ("N-IP") devices as a subset) 
responding to possibly time- varying, pulse-like input signals of finite duration, wherein 
25 the goal is to recover the integrated areas of the input signals. 

[0004] The specific embodiments described relate to processing step-like signals 
generated by detector systems in response to absorbed radiation or particles and, more 
particularly, to digitally processing such step-like signals in high resolution, high rate 
gamma ray (y-ray) spectrometers with resistive feedback preamplifiers connected to large 
30 volume germanium detectors. The method also works well with gas proportional 
coimters and scintillation detectors, to cite additional examples. The application of 
measuring the step-like output signals from y-ray detector preamplifiers to measure the y- 
rays' energies is given specific attention only because this was the area in which the 
method was first developed. 



[0005] The techniques that we have developed solve this problem generally, and 
therefore should not be construed as being limited to this specific application. Any 
detection system, for example, that produces output current signals that are integrated by 
charge sensitive preamplifiers could be treated by these techniques, whether the detected 
5 quantities are hght pulses, x-rays, nuclear particles, chemical reactions, or otherwise. The 
techniques, in fact, is not limited to "detector systems" p^r 5^, but are, in fact, general 
purpose signal processing techniques which may be broadly applied, once understood, to 
any device whose output signals that can be described in terms of multiple poles and 
zeros. The outputs from superconducting bolometers, for example, produce step-like 
10 signals that are readily treated by the invention, as do the outputs of photomultiplier 
tubes (PMTs) attached to scintillators having more than one decay time. The field of 
gamma spectroscopy, where 0.1% or less makes the difference between a bad and a good 
detector, however, provides particularly stringent tests of our techniques. 

[0006] The term "step-like signal" also requires some discussion. The output of an 
15 ideal single-pole ("1-P ") device to an ideal impulse (delta) function input, is an infinitely 
fast rise time followed by an exponential decay whose time constant Xd is characteristic of 
the pole. Viewed on a time scale short compared to Xd, this output will look like a pure 
step, while, when viewed on a time scale long compared to Xd, it will look like a pulse. A 
real 1-P device output, however, will have a finite risetime, Xr, whose duration will be 
20 determined both by the nature of the device and, particularly, by the duration of its real 
input signal. Provided that Xf is significantly shorter than Xd, a real 1-P device output 
signal, viewed on a time scale comparable to Xd, will then show a risetime region, whose 
shape may be difficult to describe mathematically, followed, after a period comparable to 
Xf, by an exponential decay with time constant Xd- The output of a device with multiple 
25 poles will be similar, except that the decay will be described by multiple time constants 
Xi, where I = 1 for the first pole, 2 for the second, etc. The output of a N-IP device will 
also be similar, with additional distortions from added minor pole or zero terms. We will 
refer to such signals, viewed on this time scale, as "step-like" or "step-like pulses." A 
pulse displaying step-like nature is shown in FIG. IB. 

30 Gamma-rav fy-rav) detection requirements 

[0007] The detection and measurement of y-ray energies is a well-established discipline 
whose primary goal is to accurately determine both the number and energies of y-rays 
emitted from some target source. The requirements of good energy resolution and high 
count rate capability usually conflict, however, since count rates are enhanced by 
35 increasing detector volume, which increases output signal distortion and so degrades 
energy resolution. High count rates also degrade energy resolution directly due to 
practical problems in preamplifier design. 
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Description of the problems 

[0008] The field of y-ray detection is highly developed. A fairly comprehensive 
introduction to the state of the art may be found in the volume "Radiation Detection and 
Measurement, 2nd Ed." by Glenn F. Knoll [KNOLL-1989]. Below we note only the 
5 issues relevant to the present invention. In the first section, we discuss how pole/zero 
cancellation errors introduce a second pole, spoiling the preamplifier's single pole 
response. In the second section, we examine how the finite input signal duration, in this 
case due to charge collection, distorts the preamplifier's output from the ideal. 

Pole/zero cancellation errors 

10 [0009] FIG. lA shows a typical solid state y-ray spectrometer comprising a 
semiconductor detector diode 7 biased by a voltage supply 8 and connected to a 
preamplifier 10 comprising an amplifier 13 with a feedback capacitor C 15 and resistor R 
17. As drawn, preamplifier 10 is a single pole circuit whose response to an impulse (delta 
function) input is A exp(-t/T2), where X2 = RC and A is the area under the impulse. 

15 Because T2 is typically of order 1 ms, which is too long for the following circuits, a 
pole/zero (P/Z) network 20 cancels the pole at I/X2 and replaces it with a pole at 1/xi, 
where Xi typically is 50 [is. Gain stage 22 then amplifies and buffers the preamplifier's 
output signal for shaping amplifier 23 which feeds multichannel analyzer (MCA) 24. 

[0010] If the time duration of the current pulse arising from the charge deposited in 
20 detector 7 by a y-ray absorption is very short compared to xi, the output of stage 22 will 
be an exponentially decaying step whose amplitude is the pulse integral and proportional 
to the deposited charge, y-ray spectrometers are therefore designed to measure these step 
amplitudes to measure the charge deposited by the absorbed y-ray. Other forms of 
radiation, including neutrons, alpha and beta particles, and x-rays behave similarly and 
25 their energies are measured the same way. 

[0011] Conrmionly, however, both the input's finite duration and the pole-zero circuit's 
imperfections distort the preamplifier's response, destroying the proportionality between 
the output step's amplitude and the deposited charge and so degrading the system's 
energy resolution. Imperfections in P/Z network 20 arise from difficulties in precisely 
30 canceling the X2 component, leaving a small residual fraction, of order 1-2%, in the output 
signal. FIG. IB shows a 5% residual X2 component for ease of viewing: an exponential 
decay signal 25 with time constant X2, input to the P/Z network 20, produces either output 
signal 27 or 29, depending upon whether the residual X2 term is positive or negative. 

[0012] These X2 residuals are particularly bothersome at high counting rates, where each 
35 signal step rides upon a X2 background from all preceding steps. As these arrive 
randomly, the resulting baseline bias also fluctuates randomly in time, which the 
spectroscopy amplifier's baseline restoration circuit cannot track well. These terms. 
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which may only be a few tenths of 1%, become a significant resolution degradation at 1 
MeV where 0.05% energy resolution is desired. 

Signal risetime fluctuations and ballistic deficit 

[0013] FIG. 2 shows a preamplifier 10 front end with a cross sectional view of the 
5 detector 7 of FIG. 1, for the common coaxial geometry. The dashed lines show electric 
field line within the detector body 30, which vary considerably with local geometry. Two 
factors cause charge collection time variations within the detector and thus risetime 
variations in the preamplifier's 10 signal output: 1) the difference between carrier 
velocities; and 2) the existence of different path lengths within the detector. 
10 RAUDORPH-1982 describes these issues. These risetime variations produce ballistic 
deficit by two paths, one direct, one indirect. The direct effect is well understood, per 
GOULDING-1988: the output filter's response varies with the time dependent shape of 
the charge arrival, being the convolution of the two. A trapezoidal filter greatly reduces 
this effect in the absence of exponential decay. 

15 [0014] The indirect effect source of ballistic deficit is due to fluctuations in charge loss 
through the feedback resistor with differing risetime, as seen in FIG. 3A with two 
risetimes, 40 and 42, where FIG. 3B enlarges their peak regions. The slower risetime 
signal loses less charge and thus is larger once charge collection is complete. Even filters 
which ignore the charge collection region are still sensitive to this lost charge effect, and 

20 relatively small errors of this size can substantially degrade resolution. For a trapezoidal 
filter, the collection time difference shown FIG. 3B produces a 0.2% amplitude 
difference (2,000 eV at 1 MeV) which will degrade ideal 1.7 keV resolution to 2.6 keV. 
Ballistic deficit errors must therefore be reduced to less than 0.05% to obtain ideal 
spectrometer resolutions at 1 MeV. 

25 Charge Trapping Losses 

[0015] Charge trapping also produces errors in y-my energy measurements since 
trapped charges are lost to the measurement. The present invention does not seek to 
address this problem. 

Scintillators with Multiple Decav Times 

30 [0016] Scintillators generally convert the energy deposited in them by a photon or 
particle into a population of excited states that then decay, emitting light as they do so. 
The total light emitted may then be taken as a measure of the energy deposited. Light 
output from a scintillator is therefore also "step-like" in the sense we are discussing. In 
addition, certain scintillators possess more than one type of excited state, so that their 

35 decay is characterized by multiple decay times. CsI(Tl), for example, has two decay 
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times of about 600 ns and 4 fis, so that its output signals formally resemble trace 27 in 
FIG. IB (the real time scale is about a factor of 100 faster, of course). 

Generalizations 

[0017] It is important to note that the pole-zero cancellation errors described above do 
5 not arise from the preamplifier's use in y-ray spectroscopy, but are a generic problem in 
low noise, charge sensitive preamplifiers. Similarly, while the described risetime 
fluctuations described above were attributed to the geometry and construction of large 
volume Ge detectors, it is clear that such problems fundamentally arise from interactions 
between the finite charge collection time and the electrical characteristics of the 

10 preamplifier and not from the physics of the collection processes. Geometric heat flow 
variations in the photon absorbing mass produce similar effects in the superconducting 
bolometers mentioned earlier. Risetime issues are therefore a potentially general problem 
as well, and may need to be corrected for in other, non-y-ray, detectors whenever the 
highest measurement accuracies are required. The methods we describe offer just that 

15 capability. Further, the terminology "single-pole" or "multiple-poles or zeros" comes 
from the LaPlace Transform treatment of differential equations describing time variant 
phenomena. Any device which shows "single-pole" behavior, for example, will display 
exponential time decay in response to an impulse input and, therefore, may be, for 
example, mechanical, thermal, chemical, or magnetic in nature in addition to the 

20 electronic case presented here. Our method can be directly applied to these devices as 
well, as will be apparent from the teachings herein. 

Existing Correction Schemes 

[0018] The prior art deals with pole/zero errors in two ways: first, by canceling X2 as 
accurately as possible; and, second, with baseline restoration schemes which try to track 
25 the shaping amplifier's "no signal" output as closely as possible, an approach which 
degrades as rates becomes high. We have not found any approaches which measure 
and/or correct for the effect directly. 

[0019] Over the years, various heuristic schemes have been developed which attempt to 
compensate for ballistic deficit. RADEKA-1982 introduced trapezoidal filtering and 

30 developed a time-variant implementation, using a gated filter following a semi-Gaussian 
shaper, that provided significant resolution improvements. WHITE-1988 proposed a 
different gated integration approach using a series switch to excise the charge collection 
region out of the preamplifier signal entirely. The final circuit was complex and had 
enhanced deadtime problems. GOULDING-1988, RAUDORF-1982, and SIMPSON- 

35 1990 disclose schemes that depend on directly measuring the signal's risetime tf and 
correcting the energy filter output by a term like tr". These approaches are complex to 
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implement and require precise expert adjustments to operate. The underlying 
assumptions are not particularly valid and improvements in energy resolution have been 
modest in practice. 

[0020] HINSHAW-1991 and KUMAZAWA-1998 describe attempts to correct for 
5 ballistic deficit by capturing peak amplitudes from two filters which respond to the 
ballistic deficit differently, one an energy measuring filter and one a differentiating (or 
bipolar shaping) filter. Typically a significant fraction of their difference in peak heights 
is added to the energy filter's peak to correct it. 

Related Art 

10 [0021] There is some related art wherein the details of the shapes of the preamplifier 
output signals are sampled digitally and used either to distinguish between different types 
of particles absorbed in the detector (e.g., MILLER-1994) or to distinguish between 
single and multiple interaction events in large germanium detectors. See, for example, 
TAKAHASHI-1994 and AALSETH-1998. 

15 [0022] WARBURTON-1997, WARBURTON-1998, and WARBURTON-1999 

describe methods for implementing digital filtering and x-ray spectroscopy. While they 
do not address the issues under consideration, some of their filtering techniques will be 
employed in the present invention and are referenced in the specification. 

SUMMARY OF THE INVENTION 

20 [0023] The present invention provides techniques for measuring a step-like output 
signal from a nominally single-pole (N-IP) device or a device with multiple poles and 
zeros in response to a pulse-like input signal to determine the integrated area of the input 
signal. The invention particularly addresses the possibility that a single-pole device 
deviates from ideality due to the presence of additional poles, zeros and/or a DC offset 

25 and that the input pulses may have finite time durations and variable amplitudes. In a 
specific example, the measurement determines collected detector charges from step-like 
preamplifier output signals in the presence of both risetime fluctuations and imprecise 
pole/zero cancellations. The invention can be implemented with digital, analog, or hybrid 
circuitry. 

30 [0024] One aspect of the invention measures the integrated area Qgj of a pulse-like 

input signal applied to a device characterized by one or more decay time constants in the 
device's impulse response function. This is accomplished by measuring a step-like 
output signal, referred to as the primary signal, provided by the device in response to the 
pulse-like input signal. In this aspect of the invention, the primary signal is directed into 

35 a plurality of signal paths to create a set of secondary signals, the secondary signal in at 
least one of the signal paths is filtered and/or delayed, a set of weighting coefficients are 



applied to the secondary signals within the signal paths, and the secondary signals are 
summed to provide a time correlated, weighted filter sum signal. In some 
implementations, the secondary signals are all summed at a single summation point, while 
in others, the secondary signals are summed at a plurality of summation points. The 
5 secondary signals are captured in at least one capturing operation after any filtering or 
delaying operations to produce a value of the time correlated, weighted filter sum signal, 
following the capturing and summing operations, that is a measure of the integrated area 
QgT of the input signal. 

[0025] Another aspect of the invention provides a method for measuring the integrated 
10 area Qgj of a pulse-like input signal applied to a device characterized by one or more 
decay time constants in the device's impulse response function by measuring a step-like 
output signal, referred to as the primary signal, provided by the device in response to the 
pulse-like input signal. In this aspect, the method comprises creating a set of secondary 
signals by directing the primary signal into a plurality of signal paths, performing a 
15 filtering or delaying operation in at least one of the signal paths, applying a set of 
weighting coefficients to the secondary signals within the signal paths, summing the 
secondary signals to provide a time correlated, weighted filter sum signal; and performing 
at least one capturing operation after any filtering or delaying operations to produce a 
value of the time correlated, weighted filter sum signal, following the capturing and 
20 summing operations, that is a measure of the integrated area Qgx of the input signal. 

[0026] Another aspect of the invention provides apparatus for measuring the integrated 
area Qgx of a pulse-like input signal applied to a device characterized by one or more 
decay time constants in the device's impulse response function by measuring a step-like 
output signal, referred to as the primary signal, provided by the device in response to the 

25 pulse-like input signal. The apparatus comprises a plurality of signal paths that receive 
the primary signal (the signals traveling along the signal paths being referred to as 
secondary signals), at least one filter and/or delay element in at least one of the signal 
paths, weighting circuitry that performs a weighting function on the secondary signals 
within the signal paths, summing circuitry that sums the secondary signals to provide a 

30 time correlated, weighted filter sum signal, and capturing circuitry that captures the 

secondary signals after the secondary signals have encountered any filter or delay element 
in the signal paths to produce a value of the time correlated, weighted filter sum signal, 
following capturing and summing, that is a measure of the integrated area Qgx of the 
input signal. 

35 [0027] In accordance with the above-mentioned aspects, the filtering, delaying and/or 
capturing operations (filter and/or delay element and/or the capturing circuitry in the case 
of apparatus) establish a defined time correlation between the secondary signals with 
respect to one another prior to the secondary signals being sunmied, and. Further the 
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weighting coefficients applied to the secondary signals are selected, based on the nature 
of any filters used in the filtering operations, on the defined time correlation established 
between the secondary signals, and on the values of the one or more decay time constants 
associated with the device, to compensate the measured area Qgj for the risetime 
5 structure of the step-like pulse, for the presence of more than one decay time constant, or 
for both. 

[0028] In accordance with the above-mentioned aspects, the order of weighting, 
filtering, and delaying may vary. For example, weighting may occur prior to filtering or 
delaying, filtering may occur prior to delaying or weighting, or delaying may occur prior 

10 to filtering or weighting. The order can differ between different secondary paths. 

Further, capturing may occur before or after summing. In some implementations, the 
capturing includes capturing one or more of the secondary signals following any filtering 
or delaying operations in their associated signal paths so that the summing of the 
secondary signals is the summing of weighted captured values of the secondary signals 

15 propagating along their associated signal paths. In other implementations, the capture 
occurs after the secondary signals have been subjected to any weighting, filtering, and/or 
delaying, and after they have been summed. 

[0029] In another aspect, the present invention contemplates processing the device's 
(e.g., preamplifier's) output signal using a set of one or more shaping filters. Where 

20 plural filters are used, they t3^ically have different time constants. A set of samples of 
the outputs of this filter set is captured in such a manner that the multiple sample values 
bear prescribed time relationships to one other. We refer to this set of sample values as a 
correlated multiple output sample (or "cMOS" for short). The relationships among the 
individual sample values may be determined by the times of their capture, by delay 

25 elements inserted in the signal paths, or some combination of the two. Further, the 

different individual sample values can be obtained from associated different filters, or 
plural sample values can be obtained from the same filter, but captured at different times 
or with different delays. The term cMOS is intended to cover these multiple possibilities 
and is discussed further in § 5.1 below. 

30 [0030] In a first set of implementations, the input pulse's integrated area is determined 
by capturing a cMOS in response to detecting a step-like signal (sometimes called an 
event). A weighted sum of the individual sample values in the cMOS (sometimes 
referred to as cMOS values) is then formed as a measure of the input pulse's integrated 
area (e.g., total charge). The weighting factors can be computed directly from information 

35 about the device's decay constants, the filter set, and the prescribed time intervals in the 
CMOS. 
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[0031] In a second set of implementations, the input pulse's integrated area is 
determined by forming a continuous weighted sum of the cMOS values and then 
capturing a single value of this composite signal as the desired measure of the input 
pulse's integrated area. In this context, the term "continuous" refers to digital and analog 
5 signals. The weighting factors are computed identically in both sets of implementations. 
The capturing operation in the second set of implementations can be either in response to 
detecting a step-like signal in the device's output (as in the first set of implementations) 
or in response to detecting a maximum in the composite signal itself. The signal paths, 
which, in the first set of implementations, contained filtering and/or delaying operations, 

10 may also contain weighting operations here. By allowing our concept of the signal paths 
leading to the summation point to include the processing of discrete (i.e. sampled) signal 
values, we see that the only conceptual difference between the first and second sets of 
implementations is whether the sampling occurs before or after the sunmiation point. In 
this generalized view, the device signal is transformed into the desired area measure by 

15 passing through a set of secondary signal paths which contain filtering, delaying, 

weighting, capturing and summing operations, and wherein the sequence in which these 
operations is performed is flexible. In the first set of implementations, the sequence is 
filtering, delaying, capturing, weighting and summing, while in the second set the 
sequence is filtering, delaying, weighting, summing and capturing. Once this concept is 

20 understood, then clearly other sequences, such as delaying, filtering, weighting, capturing 
and summing or weighting, delaying, filtering, summing and capturing are equally 
possible. 

[0032] The resultant composite filter may thus be described as being a time correlated, 
weighted filter sum-filter (TCWFS-filter), where the secondary paths produce a weighted, 

25 filtered set of correlated multiple output signals (also cMOS, since these signals in their 
sampled form are the cMOS defined above). Thus the TCWF nomenclature extends the 
cMOS nomenclature by explicitly noting the presence of the weighting and filtering 
operations in the secondary signal paths, along with the delaying operations that produce 
the time correlations. TCWFS-filtering, therefore, may either be carried out 

30 continuously, in real time, or contingently. In the continuous case (our second set of 

implementations) the weighting and summing operations operate continuously, just as the 
filtering and delaying operations do, and the desired sample is captured from the 
composite output signal. In contingent time correlated, weighted filter sum-filtering (our 
first set of implementations), only some operations, such as filtering, are carried out 

35 continuously — a set of secondary signal values is captured in response to the detection of 
an event and the remaining steps required to produce a time correlated, weighted filter 
sum-filter output value are carried out only on this captured set of values. Contingent 
operations typically include weighting and adding. However, delaying may also be 



9 



carried out implicitly by adjusting the samples' capture times. Either the continuous or 
contingent approach produces a single final value representing the input pulse's 
integrated area. Which approach is implemented depends upon the application. In 
general the continuous approach produces higher throughput, while the contingent 
5 approach uses lower cost components and is more flexible. 

[0033] The TCWFS-filter sum of the cMOS set is the composite signal that contains at 
least some signal values that are proportional to the pulse's integrated area (or the energy 
of the event that generated it), provided that the weights are properly valued in relation to 
the filter functions, the delays, and the time constants in the device. These signal values 

10 will be local maxima in the time correlated, weighted filter sum-filter signal and may be 
captured by peak sensing techniques. They also occur a fixed time interval after the 
arrival of the step-like pulse in the primary signal and so may also be captured by timing 
methods. We show how the filter time constants and other elements of the method may 
be adjusted to increase the duration of the period during which the composite output 

15 signal is proportional to the pulse's integrated area, thereby relaxing the timing 

constraints required in the time capture method. Increasing the duration also increases the 
range of input pulse risetimes that can be accommodated. 

[0034] For reasons of engineering efficiency, components may be shared when 
generating the secondary signal set. Thus, for example, the same filter, but with different 
20 delays may be used to generate two members of the set. In digitally implemented 
systems, a decimating step may precede filtering in order to reduce the resources 
required. In this case, the same decimator may precede several filters. 

[0035] Particularly in the contingent approach, the term time correlated, weighted filter 
sum is intended to cover the multiple possibilities by which appropriate time correlations 

25 may be established between the members of the captured sample set and is discussed 
further in § 5.1 below. The individual sample values captured from the time correlated, 
weighted filter set secondary signals may sometimes be referred to as the time correlated, 
weighted filter set values. The time relationships among the individual time correlated, 
weighted filter set sample values may be determined by the response times of their filters 

30 (implicit delays), by delay elements inserted in the signal paths (explicit delays), the times 
of their capture, or some combination thereof. Further, the different individual sample 
values can be obtained from associated different filters, or plural sample values can be 
obtained from the same filter, but captured at different times or with different delays. 

[0036] It is important to note that filtering, delaying, and weighting are independent 
35 operations in TCWFS -filtering, and can be carried out in any order. As a first example, it 
is this flexibility that make it easy to implement contingent TCWFS-filtering in which 
cMOS values are captured in response to a pulse and sent to a digital signal processor for 
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the weighting and summing operations. Thus, while in continuous TCWFS-filtering the 
delayed filter outputs are weighted and summed prior to the capturing step (so that the 
captured value is immediately an estimate of the input pulse's integrated area), in the 
contingent case the capturing of each individual secondary signal can be carried out 
5 independently of the others. Thus time of capture can used to introduce the proper delay 
for one filter's output while another may have a hardwired delay built in. As a second 
example, in certain systems, particularly with analog implementations, it can be effective 
to carry out weighting (i.e. gain adjustments) and/or delaying prior to filtering. Third, 
because filters incorporate implicit time delays, it is also possible to design time 
10 correlated, weighted filter sum-filters that do not require explicit time delays to achieve 
the time correlations upon which the method depends. 

[0037] For the invention method to work accurately, the used weighting factors are 
chosen to properly reflect the M-P device's decay constants, the response functions of the 
filters, and the time correlations between the different time correlated, weighted filter set 
15 signals. We show how the weighting factors can be computed mathematically, using an 
analytical description of the step-like output pulse shape and the filter time response 
functions. 

[0038] While specific embodiments of the invention use triangular, trapezoidal and 
mnning average filters, the invention does not require specific filter shapes. Further, the 

20 invention does not require that the filter's capture times be precisely located relative to 
the step-like signal's leading edge. Rather, the method derives its accuracy from 
repeatably reproducing the set of prescribed time relationships between the sample values 
in the cMOS. The underlying capture method is therefore time-based rather than 
amplitude-based, which differentiates the invention from prior art methods which capture 

25 filter samples based on their maximum amplitudes. 

[0039] While the preferred implementation uses digital signal processing, 
implementations using solely analog processing or hybrid approaches are also feasible. 
Embodiments of the invention require neither a measurement of the step-like signal's 
risetime nor more precise information about its arrival time than is found by the pileup 
30 inspection circuits of conmionly available commercial shaping amplifiers. 

[0040] When using TCWFS-filtering, however, measurement noise is increased by the 
number N of additional filter samples required to compensate for the N-IP device's 
multiple terms. For those cases where this increased noise is significant, we also show 
how N may be reduced by creating a parameterized model of the non-ideal terms; making 
35 baseline measurements to determine the parameters; and then using them to correct area 
(e.g., charge) measurements made using simpler sets of filters which have less noise. In 
some implementations, these baseline measurements are TCWFS-filter values captured at 
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times when the filters are not processing the step-like signals. In the detector- 
preamplifier case described in detail, the resultant implementation achieves the low 
electronic noise levels of conventional trapezoidal filters while also eliminating resolution 
loss due to both ballistic deficit and high count rates. 

5 [0041] A further understanding of the nature and advantages of the present invention 
may be realized by reference to the remaining portions of the specification and the 
drawings. 

BRffiF DESCRIPTION OF THE DRAWINGS 

[0042] FIG. 1 A shows a schematic of a semiconductor y-ray detector front end, 
10 including a resistive feedback integrator, pole/zero network and output gain stage; 

[0043] FIG. IB shows voltages generated by the preamplifier portion of the circuit in 
FIG. lA when the pole/zero circuit is not ideally adjusted; 

[0044] FIG. 2 expands the detector structure of FIG. lA to demonstrate the geometric 
and other factors which contribute to signal risetime variations; 

15 [0045] FIG. 3 A shows the variation in risetime between pure electron and pure hole 
collection for a simple collection path in the detector structure of FIG. 2; 

[0046] FIG. 3B shows an expanded view of the signal peaks in FIG. 3A; 

[0047] FIG. 4 shows a step-like feature in a general preamplifier output, broken into its 
component signals, and establishes a nomenclature system for describing these signals 
20 mathematically; 

[0048] FIGS. 5A and 5B show how to compute the appropriate intervals between 
capture times of the multiple filters; 

[0049] FIGS. 6A and 6B show a pair of purely digital circuits which implement the 
described technique; 

25 [0050] FIG. 7 shows a purely analog circuit which implements the described technique; 

[0051] FIG. 8 shows a hybrid analog/digital circuit which implements the described 
technique; 

[0052] FIG. 9A shows a schematic of the invention method implemented in a gate array 
as a continuously operating time correlated, weighted filter sum process; 

30 [0053] FIGS. lOA, lOB, IOC and lOD show schematics of four implementations of the 
invention method that achieve equivalent functionality through different arrangements of 
operations in their secondary signal paths; 

[0054] FIG. 10 shows input, secondary and output signals from a simple continuously 
operating time correlated, weighted filter sum filter of the type shown in FIG. 9A. 
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[0055] FIG. IIA shows generated output device signals with three different rise times; 

[0056] FIG. IIB shows the modeled response of the circuit of FIG. 9A to the signals in 
FIG. IIA; 

[0057] FIG. lie shows an enlarged view of the flat-top regions of the signals in FIG. 
5 IIB; 

[0058] FIG. 12 shows a simplified version of the circuit in FIG. 9A integrated into a 
complete spectrometer using timed peak capture; 

[0059] FIG. 13 shows a shows a different simplified version of the circuit in FIG. 9A 
integrated into a complete spectrometer using peak detect and capture; 

10 [0060] FIG. 14 shows an output trace from a CsI(Tl) scintillator/PMT detector showing 
two decay lifetimes; 

[0061] FIG. 15 shows a model of a scintillator output using two decay lifetimes; 

[0062] FIG. 16 shows the relative insensitivity of the energy resolution of a five term 
filter to count rate when working with signals of the type shown in FIG. 14; and 

15 [0063] FIG. 17 shows the normalized output traces from the detector of FIG. 14 when 
exposed to a gamma-ray, a proton, and an alpha particle of equal energies. 

DESCRIPTION OF SPECIFIC EMBODIMENTS 

1, Introduction 

[0064] Applied to preamplifiers that produce step-like output pulses in response to 
20 detector energy absorption events, the invention time correlated, weighted filter sum- 
filtering (TCWFS-filter) technique includes the following three basic steps: 1) 
generating, by applying different sets of filtering, delaying, and weighting operations to 
the preamplifier output, a set of time correlated, weighted, filtered secondary signals, or 
correlated multiple output signals (cMOS); 2) summing the secondary cMOS signals to 
25 produce a time correlated, weighted filter sum output signal; and, 3) capturing a value of 
the time correlated, weighted filter sum-filter output from which the energy of the 
detected event can be accurately recovered. The accuracy of the method relies upon 
selecting the applied weights to properly reflect the filters' transfer functions, the values 
of the time delays, and the values of the one or more decay times in the preamplifier's 
30 output response. 

[0065] Since only a single value of the TCWFS-filter output is used in the energy 
recovery step, the method can be designed to operate either continuously or contingently, 
in response to a detected step-like pulse. In the contingent approach, at least the filtering 
is carried out continuously and a set of values of the cMOS secondary signals is captured 
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in response to detecting a step-like pulse. The remaining steps to generate the TCWFS- 
filter output value are then carried out only on this set of captured values. The individual 
members in the captured set need not be captured simultaneously and, in particular, 
adjusting their capture times is explicitly presented as a means of contingently 
5 implementing necessary delays. In the continuous approach, the cMOS signals include 
weighting and are summed continuously, with a value of the summed TCWFS-filter 
signal then being captured to recover the area of the input pulse. 

[0066] The method is based upon an understanding of filters as producing a convolution 
of a filter weighting function with an input signal, so that, at any instant, the filter's 

10 output contains time correlated information about the signal over some past period. By 
capturing filter outputs at particular times (or, with delays, at the same time), we can thus 
capture information about different time intervals of the signal and then combine them to 
recover some characteristic of the signal, such as its integrated area. The appropriate time 
correlations in the TCWFS-filter method are thus achieved by adjusting the filters' 

15 response functions, inserting delays into the signal paths, setting the times of capture, or a 
combination thereof, and are independent of the specific risetime shape of any particular 
event. 

[0067] Our overall goal, therefore, is to describe how to select an appropriate set of 
filters; how to define the appropriate time delays between the cMOS filter set signals; and 

20 show how to compute the required weighting functions. Our description will move from 
the specific to the general. We will start by creating an analytic model of an isolated step- 
like pulse and describe how to manipulate summations of data values in the pulse to 
obtain the integrated area of the input pulse that produced it. From this static case we will 
move naturally to contingent time correlated, weighted filter sum-filtering, wherein we 

25 use a mixture of timing and delays to capture these same summations from step-like 
pulses in real time and then perform the same manipulations (e.g. weighting and 
summing) offline to recover the associated areas. Finally, in the last step of the 
development, we make the weighting and summing operations continuous in order to 
achieve a real time TCWFS-filter. 

30 [0068] The theory underlying the invention technique is based on modeling the device's 
(e.g., preamplifier's) step-like output response (the "output") to the input charge (or 
impulse: the "input") generated by the event as the sum of a small number of analytical 
functions which are represented by the same number of amplitudes. In a specific 
embodiment, five amplitudes are used to describe: the xi and T2 exponential components 

35 prior to the step-like signal, the ideal amplitude of the step itself, an error term due to the 
ballistic deficit, and a DC level. 
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[0069] We first reduce the number of amplitudes to four by making differential filter 
measurements which are DC level insensitive. To find these four amplitudes, we make 
four independent measurements (from three filters) in the step's vicinity, express each 
measurement as a linear combination of the four amplitudes, and solve the resultant set of 
5 four linear equations to recover the ideal step amplitude, which is proportional to the 
input charge (or impulse area). Obtaining a set of equations which is well conditioned 
and insensitive to preamplifier noise, ballistic deficit, and T2 amplitude is then primarily a 
matter of selecting an appropriate set of filters and time correlated capture points, which 
then define the cMOS. 

10 [0070] We observe that, although we are considering the preamplifier to be a 

"nominally single-pole device," we do not actually made any assumptions about the 
magnitudes of the two decay terms or of the DC component in developing the analytic 
model of the step-like pulse. The developed technique is therefore general, and can be 
directly applied to any complex pulse shape, produced by a device with various numbers 

15 of poles and zeros, for which an analytic model can be developed. 

[0071] Step amplitudes obtained using the four-filter measurement technique have 
almost a factor of 3 higher electronic noise than simple trapezoidal filters. In cases where 
this is important, we show how simpler filters can be corrected for the X2 and DC 
components by developing a parameterized model of these terms and making one or more 
20 baseline measurements to determine their values. When multiple measurements are 
made, the parameters must be compensated for their evolution in time, creating a time- 
compensated model. 

[0072] The description below is organized as follows In § 2, using a static, captured 
step-like signal, we locate several filters in the vicinity of its leading edge and 

25 mathematically describe their outputs in terms of the base amplitudes. In § 3, we use 
these results to set up and solve the set of linear equations to obtain the input charge. 
This solution has the form of a weighted sum of the filter values, showing the genesis of 
the weighting coefficients. In § 4, we evaluate the technique for two filter sets which 
compensate simultaneously for inexact pole/zero cancellation and ballistic deficit. In § 5, 

30 we introduce the concept of the cMOS to show how to implement the invention 

contingent TCWFS -filtering in real time processing situations. In § 6, we will describe 
various circuits which implement this technique: a purely digital pair, a purely analog 
circuit and a hybrid analog-digital circuit. In § 7, we analyze our technique's electronic 
noise. In § 8, we show how a simpler, three filter TCWFS-filter measurement can be 

35 used to achieve similar accuracy with lower noise, provided time recorded baseline 

measurements are made to correct for the T2 and DC preamplifier components In § 9 we 
show how continuous TCWFS-filtering can be implemented and show how it produces a 
superior flat-top filter output in the presence of multiple poles and charge collection time 
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variations. Finally, in § 10 we show how a 5 term filter produces superior high count rate 
performance with CsI(Tl) scintillators, which show 2-pole outputs. 

2. Mathematical description of generalized filter outputs 

2.1. Decay components of the step-like signal 

5 [0073] FIG. 4 shows a typical preamplifier step-like output signal as it might look if it 
were digitized by an analog to digital converter (ADC) and saved in a computer memory. 
The shown time scale is arbitrary, with zero set a little before the onset of charge 
collection for convenience. To simplify the image, only a few members of the set of 
discrete values have been identified explicitly: subscripted capital Q's denote values at 

10 specific instants since the signal really represents the charge integrated on feedback 

capacitor 15 by amplifier 13. The xi component has three regions: an exponential decay 
region 50 prior to the step (i.e., times less than 0), a charge collection region 52, and an 
exponential decay region 53 following the step (i.e., times greater than 1). The curve 55, 
showing how curve 50 would have continued to decay had the step not arrived, represents 

15 the continued time decay of all previous steps, that is, the preamplifier's first pole residual 
response to those previous steps. The signal's T2 component has the equivalent four 
regions 57, 58, 60, and 62. Finally, there is the preamplifier's DC offset value Bq 65, 
which is constant. The Xi and X2 charges injected by the y-ray pulse decay independently, 
since the preamphfier is a linear device. The total collected charge Qgx generating this 

20 event thus divides into Qigx = (l-p)QgT for xi and Q2gT = PQgT for X2, P being the 
division ratio. 

[0074] It is important to note that, in the following, we will be describing a signal 
processing technique which relies on the mathematical form of the signals and not on 
their physical source. While we develop the technique in the context of preamplifier 

25 signals from y-ray detector preamplifiers, the technique itself can be usefully applied to 
any signals which have the appropriate form. Similar outputs would generally be 
obtained from other nominally single pole devices (N-IP devices) in response to impulse 
inputs of extended duration and variable amplitude. In the present case these signals 
result from a single-pole electronic circuit with second-pole error terms responding to an 

30 input pulse from a detector, in other cases the detector itself may produce signals of the 
appropriate form. The temperature response of a superconducting bolometer, for 
example, shows this behavior [STAHLE-1999] and could be treated by the technique, as 
does the light output from CsI(Tl) scintillator crystals [SKULSKI - 2001]. 
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2.2. Running-sum filter output 

[0075] If r values are sampled at uniform intervals At from an exponential decay with 
time constant tj, then the i*'^ value Qj,i is given by: 

Qj4 = Qj,0 exp(-i At/Xj) = Qj.o bji, for (0 < i < r-1) (1) 
5 where bj = exp(-At/Tj) is the decay per At. Therefore, a digital filter, whose response 
function is defined by weighting constants { Wi,k, l<i<r-l }, acting on the values in Eqn. 1 
gives: 

ajk,r = .|^Qi,)Vi;,=QiO.|ib>i,k^QjoAjk^^ 

where the * marks this as a general filter result. This shows that the filter output is linear 
10 in the amplitude Qj,o, multiplied by A*j,k,r> which is a constant for a particular filter and 
decay constant Xi. For a simple running-sum filter (see WARBURTON-1999), w equals 
unity, so 

1-b^ 



(3) 

Eqn. 3 also defines the division factor rj, which is equal to the sum to infinity. These are 
15 finite impulse response filters whose outputs are responsive to signal inputs up to r steps 
in the past. 

2.3. General trapezoidal filter output with risetime in its gap 

[0076] Next, consider a running-sum a^H- of m samples starting at time 1 over the signal 
in FIG. 4. Decomposing the signal into 5 components with amplitudes {Qgi3» Ql3» Qg23» 
20 Q23, and Bq}, Eqns. 2 and 3 give: 

= QglsA 1 ,m + Q13A 1 ,n. + Qg23A 2,m + Q23A2,m + 0 , (4) 

where the index j in Aj^m equals 1 or 2 for Xi and X2. The equivalent sum a^., starting at 
time -4, which only sees amplitudes { Q13, Q23, and Bq} is then: 

^s- = QioAi,ni + Q2oA2^ + mBo . (5) 

25 The difference (a2+ - CJ^-) is a trapezoidal filter as of risetime m and gap g. Since Qio = 
br'n-g Qi3 and Q20 = b2-"™-S Q23: 

= ^Sf - = Qgis^ i,m + Qi3a(l.ni,.m-g) + Qg23A2.n1 + Q23«(2»ni,-m-g) ^ 

where a(j,r,s) = Aj,r(l-bjS) and the DC components cancel. In Eqn. 6, we stress that Q13 
represents, at time instant 3, the preamplifier's dominant pole (xi) response to all 
30 preceding step-like signal events; Q23 similarly represents the secondary pole (X2); while 
Qgl3 and Qg23 are their amplitude responses to the present event. Eqn. 6 takes its simple 
form because the signal rise occurs entirely within the filter gap region g, so only 
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decaying exponential terms are measured outside. Thus az is relatively insensitive to the 
precise location of the step within the gap region because its weighting constants are 
negligible (zero, in fact) there. At any time step, the output of this filter is responsive to 
signal values up to 2m+g steps into the past. 

5 2.4. Triangular filter over the charge collection region 

[0077] To explicitly probe the step's risetime region, we will use a triangular filter (gap 
= 0) Og = (ag+ - GgS) of risetime g, whose subtractive leg Og. sits precisely within the gap 
of a^. (This relationship will also be easy to achieve when we implement a digital real 
time system.) We will show that the precise location of the step within the gap is not 
10 important. Then, as above: 

= QgisA i,g + Q i3«(l,g,-g) + Qg23A2,g + Q23«(2.g^-g) -.2 q ig.i -.2 q2g4 

where the qjgj are the sampled charge values at each interval in the subtracted running- 
sum region, each comprising the new charge increment collected in that interval plus the 
exponentially decaying remnants of the charges collected in the previous intervals. 

15 [0078] Since each charge increment decays to zero independently, while their sum 
equals the total collected charge, we can write 



(8a) 
(8b) 



OO 

^2 Q2gT= ^^2gi 

Splitting the infinite summation in Eqn. 8a, however, also gives: 

20 ""i QigT= '■i QgTd - P) = iliqiga+ Qgi3 ^ 

so that 



Sqig4 = ^i{QgT(i-P)--Qgi3}-^iQci 



(10) 



The correction charge Qd defined in Eqn. 10 is the difference between the total xi 
collected charge and the amount remaining at the end of gap g. Qd will typically be of 
25 order 2% of Qgj. As a small correction term, Qd will be relatively uncorrelated to Qgi3, 
so the set of linear equations we solve for Qgi3 and Qc be well conditioned under 
inversion. 

[0079] Similarly, for the T2 component in Eqn. 8b: 
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(11) 

where Qc2 will be smaller than Qc by approximately the ratio T1/T2. Eqns. 10 and 11 
allow us to express Qg entirely in terms of Qgx, the total charge, and values from outside 
of the charge collection region, where the signals behave strictly as decaying 
5 exponentials. Equivalent integral expressions may similarly be derived in the continuous 
signal case or for filters with other weighting constants. 

2.5. Eliminating O r? 

[0080] When the gap region is much smaller than both Xi and T2, Qc2 becomes 
proportional to Qc to a high degree of accuracy. This may be seen by defining Qgi as the 
10 total Y-ray charge collected at time interval i in the collection region. Since each Qgi 
decays independently, we can write for the two components: 

,|i,(l-P)Q^br=Q,„^^ (12a) 

g-1 



.2;pQgibr=Qg23 



(12b) 



while, by definition, 



15 i = o^ " ^ ' ^^'and (13a) 



P)Qgi = (l-P)QgT. 



i=0 



Thus, for Qc: 



g-1 

Qc = ( 1 - P)QgT - Qgi3 = (1 - 3X5, Qgid - bD 



(14) 

Expanding big-> to good accuracy as (l-(g-i)At/Ti) gives (l-biS-i) equal to (g-i)At/Ti so: 

Qe = (l-P)f,|iQ,(g-i) (^3^ 

Similarly: 

Q'^=p^,li (16) 

where P2 is equal to At/T2. Comparing Eqns. IS and 16: 
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Qc2 is therefore quite small, being about 0.1% of Qc or 0.002% of Qgx- We can therefore 
express any of our filters in terms of the four quantities Qgx. Qc. Qgi3 and Qg23 using: 



.2^qig.i='^iQc 



(18a) 



5 Qgi3=a-P)QgT-Qc (18b) 

ili)^^^'^= ^2^Q<= (18c) 

Qg23=PQgT-FQc. (18d) 

3. Solving for collected charge OeT from the captured waveform 

3.1. Filter expressions in terms of the base amplitudes 

10 [0081] § 2.5 shows how all difference filters in the vicinity of the step's leading edge 
can be expressed as linear equations in the four amplitudes QgT» Qc> Qgl3 Qg23 (non- 
difference filters will also require Bq). Therefore we must make four (or five) 
independent filter measurements to solve for Qgj. While a wide variety of filters could 
be used, certain sets will be better conditioned than others and will also require less 

15 measurement precision to accommodate step risetime variations. Research and 
experimentation have shown the following set to function well: 1) the trapezoidal 
"energy" filter of risetime m from Eqn. 6; 2) the triangular "risetime probe" filter of Eqn. 
7; and, 3) a pair of triangular "xiand T2 background" filters of risetime m/2, situated at 
times -4 and +1 in FIG. 4, respectively. It is important to notice that, in this selection, 

20 only the risetime probe filter processes the step's risetime, which lies in the energy filter's 
gap and is entirely outside the range of the two other triangular filters. This placement, 
where only the risetime probe filter has non-negligible weighting constants in the risetime 
region, assures that their outputs, and hence Qgx estimate, will be insensitive to the 
precise location of the signal step within the gap, which, in real time implementations, 

25 will remove any need for great timing precision. Also, as we will show in § 9, the gap 

length must exceed the longest risetime that the filter will be expected to process correctly 
That the filters share some common starting and stopping points is primarily a matter of 
mathematical and implementational convenience, but is not fundamental to the technique. 

[0082] Substituting Eqns. 18 into Eqns. 4 through 7 gives four linear equations for ax+, 
30 as-, CTs and Og in terms of Qgj, the total new charge, Qc, the correction charge, Q13, the 
residual xi charge at time 3, and Q23, the residual T2 charge at time 3: 



20 



15 



<J2> = QgT (1-P)A i^+pA2,„] - Qc[Ai^+FA2.„j + Qi3[a 1,^ 

+ Q23[A2,iJ+mBo 

<yi^ = QgxLo] - QM + Qi3[bl'"-^A ,.^] + Q23[br"^A2^ + mBo 

^2= QgT (l-P)Ai,n»+PA2.:„] -q|a,,^+FA2^ +Qi3[a(l.m.-m-g)] 

+ Q23[«(2,m,-m-g)] 



(19a) 
(19b) 



(19c) 



'g 



■f=QgT 



(l-P)A,.g+PA2.g 



-Q, 



+ Q23 



(2-b«i)+^F(2-b|) 



0((2,g,-g) 



?13 



a(l,g,-g) 



(19d) 



We divided Og by ri so that the coefficient of Qc will approximately equal unity for 
inversion stability. 

10 [0083] Using the same methods for the outputs of three filters of length m/2 (two 

triangular filters 043 and 021 whose running-sum first points are: 043 = [time +3: time 
+1] and O21 = [time -2: time -4] and a trapezoidal filter 032 = [time +2: time -2]) gives: 

04^ =- QgT[(i-P)a(i.f .f)+Pa(2,f .f)] + Qc[a(i.?.f)+Fa(z^.f)] 
-Qi3[a(»>7.f)]-Q23[4f'?)] 

O21 =-Ql3[br~'«(l'2'2)] -Q23[b2'""^ a(2'2'2)] 
<'32 = Qgi|(l-3)A i,m/2+PA2^^ - Qc[a i.n^+FA 2.^/2] 

+ Ql3[cx(l.^.-^-g)] ^ Q23[a(2.?.-?-B)] ^^^^^ 

[0084] These filters also exclude the step's rising edge for location insensitivity. While 
this is not a requirement of our invention, it significantly simplifies both the coefficient 
20 derivations and trigger system timing requirements. 



(20a) 
(20b) 
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3.2. The four-filter solution for O^ r 
Collecting results for the four filters listed at the start of § 3.1, we have: 

tJs = QgT[(l-P)A i.m+PA2.„] - Qc[Ai^+FA2,„j 

+ Q,3[(x(l,m,-m-g)] +Q23[a(2,m,-m-g)] 

■(l-P)A,,g+3A2.g 



(21a) 



-^=QgT 



-Q< 



+ Q 



13 



a(l,g,-g) 



+ Q23 



(2-b^)+^(2-b|) 
^1 



a(2,g,-g) 



(21b) 



^43= -QgT 



021 =-Q 



i-P)a(i.f .f)+pa(2.f .f)] + Qc[a(i.?.f)+Fa(2.^.f)] 
-Q,3[a(i.f.^)]-Q23[a(2.f.l)] 
br^s a(i,^ m| _ Q23[b-2'"^ 0(2,^ 



(21c) 
(21d) 



as four linear equations in four unknowns. The unknown parameters are of three types: 

1) (type Pi) the unknown charge QgT (the area of the input impulse in the general case); 

2) (type P2) charges describing the step's finite arrival time (Qc) and the residual charges 
in the Ti and T2 decay modes (Q13 and Q23) from previous events; and 3) (type P3) a set 
of fixed parameters {xi, T2, P, m and g} which describe the preamplifier's (or other 
device's) transfer function and the applied filters. Since the coefficients in [ ]'s are 
functions only of the type P3 parameters, they are constants which need be computed only 
once for a given preamplifier and filter set. 

[0085] Eqns. 21 can be solved by several methods, including casting them into matrix 
notation and then inverting the coefficient matrix. Thus: 
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_g 

^21 



(l-P)Ai,„+PA2., 
(l-P)A,,g+PA2,g 



_(2-b^)-^(2-bi) 



a(lan,-m-g) 



«(l.g.-g) 
I"! 



a(2jn,-in-g) 

a(2^.-g) 
r. 



-(1-P)a(i.f .f)-Pa(2,^.f) a(i.f.f)+Fa(2.^.f) -a(i,f .f) -"(^.f -f) 

0 0 -b7-^a(i.f.f)-b-'»-%(2.f.f) 



QgT 

Qc 
Ql3 

Q23 



becomes 
so that 

and 



QgT = Jll J 1277 + Jl3 <'43 + Jl4 ^21 . 



a=J.Q 



Q = J-*«a 



T-l 



■1(22) 
(23a) 
(23b) 

(23c) 



[0086] As a solution for Qgx, Eqn. 23c has the advertised form: a weighted sum of the 
four filter outputs. The weighting coefficients Jy-l are computed only once for a 
particular set of filters and P3 parameters {xi, T2, P, m and g} and then multiplied by the 
four filter outputs as, <7g/ri, G43 and 021 to obtain the total pulse charge Qgx (and hence 
10 the energy of the absorbed Y-ray) for each captured event. We note that if more filters are 
measured than there are parameters Q, the matrix can be inverted using least squares 
methods to obtain Eqn. 23b. 

4. Evaluation for two filter sets 

4.1. A four-filter set with typical parameters 

15 [0087] We now evaluate a typical case of the filter set leading to Eqns. 22 and 23: we 
set m equal to 4 (is, g to 1 |lis and assume typical values of Xi = 50 |is and X2 = 1 ms with 
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10 



15 



20 



25 



30 



a division ratio P = 0.04. In this case, for 40 MHz sampling, m = 160, g = 40 and we can 
compute 



J = 



SO 



and 



1.54041E+2 -1.54138E+2 -8.08795E+1 -8.00411E+1 
1.98087&-2 --1.06150E+0 -2.00007E-3 -1.99953E-3 
-2.95906EH-0 3.07604E+0 -1.53785E+1 -1.59682E+1 
O.OOOOOE+0 O.OOOOOE+0 -1.69959E+1 -1.60483E+1 



6.309486^3 -9.623578E^1 -1.593231E-2 ^1.549595E^2 
1.190638E-^ -9.600143E-1 -2.284270E^ -2.469335E^ 
1.194301^2 6.874687E^2 6.221807E-1 -6.786521E-1 
-1.264822^2 -7.280624B-2 -6.589193E-1 6.564133E-1 



(24a) 



(24b) 



Qg^, = 6.309486E^3 a2-4.811789&^ - 1.59323 1E^2 a^^ - L549595E^2 a^^ ^24c) 

Cg is typically of order 0.25a2:> so the Gg term is typically less than 2% of Qgx- In spite 
of its small size, however, Gg carries the burden of the ballistic deficit correction, since it 
is the only filter that probes the step's risetime directly. 

[0088] The coefficients in Eqn. 24c also depend upon the value p. When the pole-zero 
is perfectly balanced, (P = 0) Eqn. 24c becomes: 



QgT = 6.309486&-3 a^^- 4.811789Er^ - 1.59323 lE-2 - L549595Er-2 



(25) 



where Og's coefficient has changed by 4%, a^'s by 0.04%, and the 043*8 and aai's by 
about 0.008%. The need for P accuracy is therefore a function of how accurately one 
wishes to correct for risetime fluctuations and pole-zero inaccuracies. The J-1 weighting 
coefficients are found to vary linearly with p, which allows P to be determined 
experimentally if it is not precisely known. To do this, first select a set of P values 
bracketing its expected value and compute Eqns. 24 to produce a set of weights ^^1° 
(the coefficients in Eqn. 24C) for each P value. Next collect data for each weight set 

and measure the standard deviation aQgT(P) in Qgx (i.e., the energy resolution) 
which is obtained. Then, by plotting or interpolating the oQgT(P) results versus P, select 
a P value which minimizes aQgT(P). If the preamplifier is stable, P should only need to 
be determined once. 

[0089] We tested the accuracy of Eqn. 24c in the presence of both risetime variations 
and a sizable T2 decay term (equal to the Xi step amplitude) using a spreadsheet to 
compute the filter responses for modeled waveforms. Table 1 presents these results. As 
may be seen, the Eqn. 24c is insensitive to both risetime fluctuations and X2 components 
at the few ppm level and is three orders of magnitude better than a simple trapezoidal 
filter. 

Table 1: Modeled charge Qgx for 5 filters in the presence of varying risetimes and T2 
amplitudes for the case of P equals 0.02. 
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Risetime = 450 ns Risetime = 950 ns 



Model 


<T2> = 0 


<X2> = 1 


<T2> = 0 


<T2>= 1 


Triangular 


0.861966 


0.859787 


0.831149 


0.827170 




(100.00%) 


(99.75%) 


(96.43%) 


(95.96%) 


Trapezoidal 


0.867397 


0.862421 


0.863901 


0.858928 




(100.00%) 


(99.43%) 


(99.59%) 


(99.02%) 


Section 4.1 


0.999997 


0.999997 


1.000001 


1.000001 


(Eqn. 24c) 


(100.000%) 


(100.000%) 


(100.000%) 


(100.000%) 


Section 4.2 


1.000119 


1.000113 


1.000028 


1.000022 


(Eqn. 33) 


(100.000%) 


(99.999%) 


(99.991%) 


(99.990%) 



4.2. A three-filter simplification which still has good performance 

[0090] For cases which do not demand the precision of Eqn. 24c, or have small T2 
components, we show a simpler, three-filter set which still has good performance. We 
5 use trapezoidal/triangular filters for their insensitivity to the "DC-like" slow variation in 
the T2 component: namely the same "energy" and "risetime probe" filters as in § 4.1 and 
a second trapezoidal filter 032 (Eqn. 20c), aligned to exclude the risetime region, to 
estimate Q13. Thus, assuming P and Q23 are all zero: 























-a(l,m,-m-g) 


QgT 






0 


1 


^_fia(l,g,-g) 


Qc 


Wa32 




WA, m 


-WA, m 


wa(i,f.^-g) 


Ql3 



(26) 



10 In Eqn. 26 we reversed the sign of as in order to create a filter that is differential on both 
sides of the step's rising edge and we weighted 032 by W to determine the value that 
gives the best results. We find the lowest noise is for W = 1, so we can repeat Eqns. 24a 
and 24b to obtain: 

QgT, = 9.338925E^3 - 5.000000Er^ + 3.106015E-2 ^27) 

15 [0091] Table 1 shows that this result is sensitive to the T2 component at the few ppm 
level and to risetime variations at the 100 ppm level. For many practical cases, this is 
sufficiently accurate. Assuming a random distribution of risetimes, for example, this 
would add less than 100 eV in quadrature to the energy resolution at 1 MeV, which would 
be scarcely perceptible. 
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5. Correlated multiple output sampling (cMOS) in contingent TCWFS-filtering 

5.L Defining '^correlated multiple output signals - cMOS" 

[0092] The filter weighting Eqns. 25 and 27 derivations were developed as if the 
signals had been captured into memory, so the running-sums could simply be taken over 
5 the desired sample values. This implementation is feasible at lower data rates, but in a 
preferred implementation, as taught by WARBURTON-1997, a gate array is used to 
realize the difference of running-sum filters, which run in real time, each outputting a new 
value a fixed number of clock cycles after a new ADC sample is supplied. In general, 
then, there may be no single instant when a set of filter signal values can be captured 

10 which will reproduce the inter-filter time relationships required by Eqns. 25 and 27. 
Therefore, the outputs of the different filters will either have to be captured at different 
times, and these times will have to be carefully correlated if Eqns. 25 and 27 are to 
produce accurate values of Qgx, or appropriate delays will have to be applied so that all 
the values can be captured simultaneously. Because multiple filter signal values must be 

15 captured, and their times of capture must be properly time correlated, we call the resultant 
set of captured values a "correlated multiple output sample set," or "cMOS," and refer to 
these values as "cMOS values." Since the cMOS values are a snapshot of the delayed 
filter output signals (the secondary signals) we will also refer to these "correlated multiple 
output signals" as the "cMOS" as well. Applying the weights of Eqns. 25 and 27, we 

20 obtain a set of "time correlated, weighted filter" values. 

[0093] We also observe that the number of cMOS values required by Eqns. 25 and 27 
and number of realized filters needed to produce them need not be the same. For 
example, in Eqn. 25, the 021 and 043 filter values both come from triangular filters of 
length m/2 with their starting points separated by the interval m+g. Thus two outputs 
25 from the same realized filter can be captured to supply the two required values. 

[0094] While we have noted that the filter ranges do not have to be precisely positioned 
relative to the step 's rising edge, the accuracy of their positions in time relative to one 
other in the cMOS capture determines the accuracy that the technique can attain. This is 
because the running-sum ranges, which translate directly into filter lengths and time 

30 intervals between starting locations, are encoded into the J-^ weighting coefficients 

through terms like a(j,r,s) = Aj,r(l-bjS) (see, e.g., Eqn. 6), where r is the filter length and 
s the interval between the sum start points. For simplicity, we aligned the various filters 
so their starting points coincided, but other cases are trivially treated by scaling Eqns. 4-7 
by terms like biS where t is the offset between starting locations. Therefore, there is an 

35 intimate relationship between the weighting coefficients and the filter ranges which, if 
violated, will introduce inaccuracies. The J-*'s are computed for a particular set of filter 
sum ranges. If another set of ranges is captured, then expected cancellations between 
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filter terms will not occur exactly and the differences will appear as errors in the detected 
charge. In digital implementations, the fact that the signal samples are discrete and 
enumerated makes is relatively easy to maintain accurate time relationships by measuring 
time in "ticks" of the digitization clock and then simply counting out the required time 
5 intervals. 

[0095] The requirement of appropriate time correlations between the cMOS values 
clearly differentiates the present method from prior art as represented by the approaches 
of HINSHAW-1991 and KUMAZAWA-1988. While two filter outputs were captured 
and weighted in those methods, the capturing steps used peak detecting means. The 
10 methods are therefore amplitude based rather than time based, as in the present method, 
and the captured values carry different information about the step shapes. Further, the 
method of capturing based on peak detection assures that there will be no fixed time 
relationship between the captured samples since the two filters' times of maximum 
amplitude depend differently on step shape, which will vary from step to step. 

15 [0096] In light of the above, it is clear that there will be a wide variety of possible 

approaches to capturing filter cMOS values to form the TCWFS-filter with the prescribed 
time relationships in response to detecting a pulse-like feature in the signal. The outputs 
of single filters can be captured multiple times. Filter outputs can be captured as they 
appear from the gate array or can be delayed for convenience. We will show several 

20 examples in § 6. Which method is best will depend upon implementation details and will 
vary from case to case. We therefore intend the phrases "cMOS" and ^TCWFS-filter'' to 
cover the multiplicity of possibilities for capturing, in a real-time implementation, filter 
values whose time relationships are correct for the set of J'^ weighting coefficients that 
will be applied to them to recover the input charge from the detected event. When the 

25 capture is made in response to the detection of a step in the preamplifier's output signal, 
with the intention of using it to determine the charge deposited in the detector (in the 
general case the impulse area) and hence the gamma-ray's energy, we will often refer to it 
as an "energy (or "charge or area") cMOS or TCWFS capture" (as appropriate). 

5.2. Timing considerations 

30 [0097] Filters, in real time implementations, require finite times to compute. For 

example, consider a trapezoidal filter of length (risetime) m and gap g over a set of values 
{qi}, where qn is the last sample value ("end sample") to enter the filter: 

n n-m-g 

a^Cn) = S qi- 2 

' i=n-m i=n-2m-g ^ (28) 

If implemented in gate arrays according to the method of WARBURTON-1997, the 
35 output a2:(n) will appear 3 filter clock cycles after sample qn is presented. This is its 
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propagation delay. For a 4 fis peaking time (filter clock equals 200 ns), the propagation 
delay is 600 ns. These delays must be included in the design to assure that the cMOS 
capture constraints are satisfied. 

[0098] FIGS. 5A and 5B show how to compute sampling times for the Eqn. 24c filter 
5 set. The calculation has two steps: finding the time delays between their end samples of 
the different filters and then correcting for their different propagation delays. For the first 
step we use FIG. 5A, which shows the filter regions relative to the preamplifier output 
step 72. ax 74, C43 76, 021 77, and Og 79 were all described earlier. Ot 81 is a fast 
channel timing filter, as described by WARBURTON-1997, which is used to detect 

10 step-like signals in the preamplifier's output. The timing filter's output triggers a digital 
discriminator when it is located approximately as shown with respect to the leading edge 
of the signal step 72. We can now measure end sample delays relative to at 81 as shown, 
with the filters a^ 74, 043 76, and Og 79 appearing at times At^, At43, and Ats after at 81 
and 021 77 appearing At2i earlier. From FIG. 5A, if at has a peaking time of T, then Atx 

15 equals (m+g-2T). 

[0099] It is important to note that, while the firing of at may jitter somewhat depending 
upon the step's risetime shape, this only causes the step location to jitter within the gap g. 
The relative timings of all the filters to at, and thus to each other, are fixed by 
construction. Thus, provided Og is wide enough to catch the step for any reasonable 
20 discriminator triggering, the filter will work as designed. This insensitivity to the timing 
accuracy with which the step-like signals are detected is a significant advantage of the 
present invention. 

[0100] FIG. SB shows the second step: correcting for propagation delays. We note 
that these delays may exceed a microsecond in long time constant filters and typically 
25 increase with filter length. FIG. 5B shows propagation delays with respect to the ADC 
sample time tg 83 for a 40 MHz implementation with g = 1 |iis and m = 4 ^is, so 5tt = 75 
ns, 5tg = 5t2i = 5t43 = 300 ns, and St^ = 1,200 ns. Combining FIGS. 5A and 5B, we 
should sample the output of a given filter ax, at 

tx = ^t^+St^-Stt (29) 

30 following the discriminator firing on at's output to assure the cMOS temporal 

relationships required to extract Qgx from the captured TCWFS values. If tx is positive, 
then a counter can be used to measure tx and capture the filter value. If tx is negative, 
however, then its propagation delay 5tx will have to be increased by at least |tx| to meet 
the demands of causality (capturing must follow step detection in time). In digital 

35 implementations this is a simple operation using a FIFO memory. In our present 

example, no additional delay is required since a typical timing filter peaking time is 100 
ns, so At2i = -200 ns and 
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t2l = At2i + 5t2i + 6tt = -200 + 300 - 75 = 25 ns, (30) 

which is positive. 

6. Circuits implementing the cMOS contingent TCWFS-filtering method 

6.1, Fully digital filter implementations 

5 [0101] Two general approaches are possible for capturing the cMOS values following 
(i.e., contingent upon) step detection: 1) asynchronously capturing values as they appear 
from the sub-filters, per the discussion of Sect. 5.2, and 2) extending the filters' 
propagation delays so that they can be synchronously captured. In the latter case all 
pileup inspection is carried out before the capture signal is given, while in the former case 
10 the values are captured as they appear and the pileup inspection is used to decide whether 
anything should be done with them afterwards. We will demonstrate one implementation 
of each to show the principles. 

6.1.1. Asynchronous capture implementation 
[0102] FIG. 6A shows a preferred implementation of Eqn. 24c. The digital processor 

15 90 comprises an ADC 92, a real-time digital processing unit (RT-DPU) 93, and a digital 
signal processor (DSP) 95. Per the method taught by WARBURTON-1997, the RT-DPU 
is implemented in combinatorial logic (e.g., a gate array) and processes data at the ADC 
output rate while the DSP processes data at the event rate (impulse signal rate). The RT- 
DPU produces values a^, 043, 021, and Gg, while the DSP multiplies them by the 

20 coefficients J-^ to obtain Qgj. 

[0103] The RT-DPU 93 also has two parts, a fast channel 97 and a number of slow 
channels 98, The fast channel contains the fast timing filter at 100, a digital discriminator 
102, some pileup test and timer logic 103, a clock 104, and an output register 105. Each 
slow channel comprises a triangular or trapezoidal digital filter and an output register 

25 with either a digital delay line (if Eqn. 29 is negative) or a timer (if Eqn. 29 is positive). 
Notice that 043 and 021, both being of length m/2, are being implemented using a single 
triangular filter 04321 107. 

[0104] In operation, at 100 processes ADC 92 data until a step-like signal arrives. 
Discriminator 102 detects this step and signals the pileup inspector 103, which strobes a 

30 trigger line 106 and initiates pileup inspection. If the step's arrival time is important, the 
output of clock 104 can be captured to output register 105 by the same strobe signal 106. 
Filter 04321 's 107 output, delayed by t2i using 21 -delay 108 if required by Eqn. 29 
causality, is captured immediately in the 21 -register 109 by this trigger strobe, which also 
starts the 43-timer 110, Z-timer 111, and g-timer 112 to time 143, and tg. Since tg is 

35 shortest, the g-timer 112 times out first, capturing filter Og's 115 output in the g-register 

29 



116. Similarly, the 43 -timer 110 and Z-timer 111 time out at t43 and ts to capture filter 
<y4321 107 and as 117 outputs in the 43-register 118 and S-register 119. If this step is not 
piled up, then all four filter values are ready to be read into the DSP 95 via data bus 121. 
The pileup inspector can be implemented as taught by WARBURTON-I999 relative to 
5 the length of the filter and will not be discussed further. 

[0105] As per the cMOS discussion in § 5,1, the invention method does not specifically 
require that the number of filters equal the number of filter values to be captured. Here 
we have captured G4221 twice. 

6.1.2. Svnchronous capture implementation 

10 [0106] FIG. 6B shows an implementation of Eqn. 27 wherein the three filter values are 
captured just as the value emerges from the slowest (a^) filter. The primary physical 
difference between this implementation and that of FIG. 6A is that the timers 108, 109, 
and 110 have been eliminated and a second delay (g-delay 125) has been added to 
compensate for the differences between tz and tg and 132 from Eqn. 29. All the 

15 remaining parts are the same and carry the same reference numbers as is FIG. 6A. Notice 
that the trigger line 106 now connects directly to the output registers, rather than through 
timers as before. For this topology, the functional operation of the pileup test and timer 
103 is essentially identical to that presented by WARBURTON-1999. 

[0107] A disadvantage of this circuit, compared to FIG. 6A, is that relatively long 
20 delays may be required, which may consume excessive gate array resources. Recent 
programmable gates, however, are appearing with blocks of dedicated memory, so this 
restriction is becoming less significant. One significant advantage of this design is that, 
because of the synchronous capture, the pileup inspection logic design is much simpler. 

6.2. Fullv analog implementation 

25 [0108] While our preferred implementation is fully digital, the method is not inherently 
digital per se. In certain cases, particularly for very fast step risetimes and very short 
filter time constants, or to allow a high degree of integration, it may become preferable to 
use an analog implementation. FIG. 7 shows an analog version of FIG. 6A. The 
topologies are identical up to the point where sample and hold circuits replace digital 

30 registers to capture the filter values, so the operation of the circuit up to this point will be 
clear to one skilled in the art, given the teaching presented herein. The derivation of the J 
matrix, of course, will have to proceed using weighting functions appropriate for the 
analog filters used. Then analog multipliers (160, etc.) and an analog adder 185 are used 
to implement the J-^ coefficients and addition of Eqn. 24C. We note that, for an 

35 application with fixed filters, the multiply-and-add function is just an op-amp circuit with 
appropriate resistor inputs to its summing node. If the pileup tester 143 deems the value 
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good, the gate and shape circuit 145 passes out the adder 185 output in a form suitable for 
multichannel analysis. 



6.3. Hybrid implementations 

[0109] Shaping filter characteristics which are difficult to implement digitally may be 
5 required in some situations where the other benefits of digital processing are still desired. 
In these situations hybrid implementations may be considered. FIG, 8 shows one such 
implementation 190. The topology is nearly identical to that of FIG. 6A. The major 
differences are: 1) that the fast channel 191 is implemented in analog circuitry identical 
to that of FIG. 7 (although the digital circuit of FIG. 6A could as easily be used); and 2) 
10 that, in the slow channel section 192, the digital filters 04321 107, as 117, and Og 115 
have been replaced by analog filters 04321 1^3, aj; 194, and Og 195 followed by ADCs 
196, 197, and 198. The time behavior of this circuit is the same as was described for the 
asynchronous digital circuit in FIG. 6A. 

7. Noise considerations 
15 [0110] The Eqns. 24C and 27 filters have higher noise than simple trapezoidal filters 
since they weight the sampled set of data points with values larger than unity. For 
example, by breaking the 2m+g filtered region of Eqn. 27 into its 4 major sub-regions of 
length m/2 (ignoring the gap g for this argument) and calling running-sums over these 
regions ai through 04, we can recompose Eqn. 27 as: 

^- 9.338925E--3 + 2. 172123E-2 - 2.172123E-2 G2 + 9,338925E^3 ^^i) 

If the noise standard deviation per sample is a, then the total noise sig(QgT) will be 

,0.5 



siglpgyj = ^160 ^8.722E^ 5 a^j + 160 ^4.718Er^ a^jj = 0.299 a 

which can be compared to a trapezoidal filter of the same length: 

,0.5 



(32) 



-g(Q3=[l6o(iio)%l^^ =0.112a 



(33) 

25 Thus the Eqn. 27 filter will have approximately 2.7 times as much electronic noise as a 
simple trapezoidal filter. The result for the Eqn. 24c filter is similar. For a 4 |xs filter, 
whose electronic noise can be as low as 150 eV, this effect is significant at energies like 
100 keV (Fano noise 415 eV) but not at 1 MeV (Fano noise 1,300 eV). And, as peaking 
times shorten to obtain higher count rates, the problem will worsen at higher energies. 

30 We therefore consider how to obtain the benefits of correcting for both ballistic deficit 
and pole-zero errors without paying the additional filter noise price. 
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8. Treating non-ideal terms with baseline measurements 

[0111] Our approach is to extend the baseline averaging concept of WARBURTON- 
1997 to treat time varying terms as well. In the particular case at hand, we would like a 
way to measure Q23 using baseline measurements in spite of the fact that it both exhibits 
5 exponential decay and is increased each time an event occurs. 

8.1. Derivation of an example filter 

[0112] We begin by developing a simple cMOS TCWFS- filter, composed of three 
running averages, of lengths m, g, and m, which we make contiguous for the practical 
reasons of simplifying the math and minimizing deadtime, but not for any fundamental 
10 reason. As in previous sections, the length m filters will primarily serve to determine 
QgT, while the length g filter will be used to eliminate ballistic deficit effects by probing 
the signal step's risetime region. Starting with Eqns. 19a, 19b, and 19d, these filters are: 

<yi^ = QgTAm- QcAm + Q I3A i,m + Q23A2,m + ^^0 ^ (34^) 

<y^= +Ql3b";^Aig +Q23b2^A2,g +gBo^ ^34^^ 

15 and ""^^ +Qi3br^Ai^ +Q23b-"^A2^ ^mB^^ ^^^^ 



where 

Am = (1 - P)A i,m + PA 2,m, A^ = A + FA2,m, f = ri + Frj 



(34d) 



[0113] In Eqns. 34 we now have four types Pi of parameters: types 1 and 2 (Pi, P2), as 
before, which include Qgx, Qc, and Q13; type 3 (P3), Q23 and Bq, which describe the 

20 preamplifier's "non-ideal" behavior, Q23 being the amplitude of its second pole's 
response to any previous steps and Bq being the DC offset; and type 4 (P4): all the 
remaining fixed parameters (e.g., Ti, T2, P, m, g) which describe both the filters' and the 
preamplifier's transfer functions. We will treat Qgj, Qc> and Q13 as unknowns to be 
recovered from a particular TCWFS-filter measurement, with the intention of 

25 independently determining Q23 and Bq from baseline measurements. Eqns. 34 are simple 
enough to be solved by direct elimination of unknowns to get: 



QgT = A-^ 



where 



gj -A-'[lQ23 + MBo] 
J a 



(35a) 



1 r A . 1 



i,m , (35b) 
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L = A 



24n 




m+g 



and 



M= m 



l,in 



.m+g 



(35c) 



m 



l.g 



1^ 



(35d) 



The constants K, L, and M depend only on type 4 parameters, which are assumed to be 
time independent for a particular detector-preamplifier pair and set of filter lengths. Eqn. 
35A is therefore a transform function (the first bracketed term) between the Eqn. 34 
TCWFS detector charge Qgx, minus an error term t^Qer + MBq] ^ ^ ^j^^^ -g ^j^^ 
transform's response to the non-ideal (NI) preamplifier terms Q23 and Bq. Using Eqn. 
35a thus requires that Q23 and Bq to be determined at a time appropriately correlated to 
the cMOS capture (e.g., time 1 in FIG. 4) and substituted into the error term to produce 
an accurate value of Qgx- The cMOS capture itself can be physically implemented using 
any of the approaches described above. The subscript q on the bracketed weighted sum 
of CMOS values indicates that this was an "area" cMOS TCWFS-filter capture of the area 
filter set to determine deposited charge. 



8.2. Single measurement error term correction 

[0114] A simple, although not the lowest noise, method to estimate the error term is to 
measure it locally for each cMOS energy capture by making the same Eqn. 34 cMOS 
filter capture in close time proximity but without a step-like signal present in the filters, 
so that QgT is zero. Then, from Eqn. 35a, with Qgj equal to zero, we have: 



m 



LQ23 + MBq 



Jb 



m 



A 



g- 



, (36) 

where the subscript b is for "baseline," called so since, as shown, this "baseline" cMOS 
filter capture produces the transform's baseline response to the preamplifier's non-ideal 
terms. We may also call such a measurement a "null charge" or "null energy" ("null 
area" in the general case) capture since it is made when QgT is zero. Substituting Eqn. 
34c into Eqn. 33a gives: 



QgT=A-* 



A. 



A 



r g- 



(37) 



where the subscript q labels the TCWFS-filter charge measurement to be corrected. 
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[0115] This correction is moderately simple and has been found to work acceptably in 
some cases. However, if Q23 is large, as may be true at high output step rates, then time 
decay effects in [LQ23 + MBo]b may still be significant. This can be approximately 
corrected by measuring the time interval 5t = te - tb (in clock ticks At) between the energy 
and baseline captures and correcting by the decay term b2^^ to obtain: 



QgT = A 



T -1 



m 



As 

A 



-02 



A. 

A 



^2+ 



(38) 



In this method, if Bq is truly a constant and the time interval 5t between the energy and 
baseline captures is fixed, then the error term MBob2^'^, will also be a constant and merely 
produce an offset in the spectrum. 

10 [0116] Eqn. 38 has two significant problems: 1) increased noise, due to the noise in the 
single baseline measurement, which is the same size as in the measurement itself and 
degrades resolution by sqrt(2); and 2) it increases deadtime by at least a factor of two if 
the two measurements are independent and so separated by at least 2m+g in time. It is 
therefore useful to be able to make multiple, independent measurements of Bq and Q23 

15 and average them in order to reduce their variance. We present a method for doing so 
below. 



8.3. The time compensated baseline model 

[0117] In making filter derivations throughout this work, we have assumed that we have 
an accurate preamplifier model in terms of its several decay constants and DC offset. We 

20 therefore take the logical step and use precisely the same model to correlate 

measurements made at different times so they may be averaged. Thus, if we measure 
Ql3, Q23, or Bo at one time, then their values At clock cycles later will be given by 
bi^K5i3, b2^^Q23, and Bq, respectively, unless a step arrives within the interval At. In that 
case, at the end of the step arrival time (e.g., time 1 in FIG. 4), we increment Q13 by (1- 

25 P)QgT» Q23 by PQgT, and Bq by zero, after which their exponential decays will continue 
as before. For the purposes of baseline averaging, then, Q13, Q23» or Bq become "time- 
compensated" model (or *TC-model") parameters, whose averaged values are <Qi3>> 
<Q23>» and <Bo>, respectively. 

[0118] If we make measurements of Bq and Q23 from time to time, we can now update 
30 our estimate of Bq simply by using the method of WARBURTON-1997 or some other 
averaging scheme. For Q23 our TC-model now gives: 



Q23(i) = bt'Q23(i-l) 



(39) 
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so, if <Q23>i-i represents an "averaged" estimate of Q23(i-1), last measured at time ti.i, 
we update <Q23>i> using a new baseline measurement of Q23(i) at time ti, by the 
following extension of Warburton and Hubbard's Eqn. 7, where Ati = tj - tj-i: 



Q23(i) 

N 



(40) 



Equivalently, we update <Q23> for the arrival of an impulse signal of charge QgT(i) at 
time ti using: 

<Q23>i= bt'><Q23)(^l) + PQgT(i) (41) 

If required, Q13 can be similarly extended from its last time of measurement by replacing 
b2 with bi in Eqn, 39. 

[0119] <Q23> and <Bo> are used with an charge cMOS TCWFS-filter capture to 
determine QgT using Eqn. 35a as follows. First, <Bo> simply replaces Bq in the 
equation. <Q23>, however, has to be extended from its time of last update ti-i to the 
present capture time ti using Eqn. 39, according to our TC-model, where Ati = ti - ti-i: 



QgT = A-* 



m 



^r'^g- 



(42) 



Because <Q23> and <Bo> can be determined by averaging over many baseline 
measurements, the noise contribution of the correction term in Eqn. 42 can be made 
much smaller than in the single baseline correction procedure of § 8.2. Since the 
weighting function of ag_ is small compared to unity, the measurement noise in Eqn. 42 
will then be determined primarily by a2+ and 05:., and will be approximately the same as 
for a trapezoidal filter of the same peaking time and about a factor of sqrt(2) better than 
by using Eqn. 38. Furthermore, because a baseline cMOS or TCWFS-filter capture can 
now be made at any time the interval between successive steps is adequately long, rather 
than requiring one in conjunction with each charge cMOS filter capture, the method's 
deadtime is cut in half: to that of a single cMOS filter capture. 
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8.4. Measuring the non-ideal terms Qj^ and B o 

[0120] An effective way to measure Q23 at time ti is by making a baseline cMOS 
TCWFS-filter capture, per Eqn. 36, with the difference now that the times ti need no 
longer be correlated with charge cMOS captures. Using an averaged value <Bo> for Bq, 
plus computed values for L and M, we can then solve [LQ23 + MBoJb for Q23- 



Q23 =f 



A. 

A. 



r g- 



Jb 



(43) 
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and use it to update <Q23>i using Eqn. 40. 

[0121] We normally expect Bq to be stable and nominally rate independent, so that 
<Bo> can be measured with no input counts to the detector and used thereafter to obtain 
Q23 values. However, to compensate for any temperature and/or rate related changes in 
Bo, it is prudent to also measure Bq from time to time and update our estimate of <Bo>. 
One approach is to use pairs of baseline measurements, separated in time by At, where, if 
we denote the error term E as: 



Bo,i measurements can be averaged, as noted above, to obtain <Bo>, a low variance 
estimation of Bq, for use in Eqn. 35a. From Eqn. 45, it is clear that the longer the time 
interval At, the more accurate the measurement of Bq will be, since both the numerator 
and denominator are small differences of large numbers. In any case, enough values Bo,i 
can be averaged so that errors in <Bo> do not contribute substantially to the error in Qgx- 

8.5. Dealing with model parameter drifts 

[0122] For this compensation scheme to work well, the P4 parameter P needs to be 
precisely known and stable over time. Else both the extrapolations of Q13 and increments 
in Q23 will incorporate systematic errors which will propagate into the Qgx calculations. 
The proposed Q23 baseline measurements offer a method to monitor P and either refine 
its value or track it in time. 

[0123] In Eqn. 40 we update our average <Q23>i = b2^^<Q23>i-l at time ti with the 
measured value Q23(i)- However, we can also record the differences between these two 
values and compute their mean <AQ23> over time. If <AQ23> is not zero, then a 
systematic error exists and we can adjust p using: 



where the slope Sp = d<AQ23>/dP is previously obtained from calibration runs, made at 
system setup time, wherein P is systematically adjusted and the resultant <AQ23> values 




(44) 



(45) 



5P = <AQ23 >/Sp, 



(46) 



recorded, so that Sp can be calculated. Similar methods can be used to monitor and 
correct other P4 parameters such as Xi and T2, if appropriate. 
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9, Circuits implementing the continuous TCWS-filtering method 

[0124] The contingent method of capturing a weighted signal sum value presented 
above has been tested and applied in a variety of situation and found to work well. 
However, the time required to pass the multiple secondary signal samples to the digital 
5 signal processor, plus the time required to form the weighted signal sum, has been found 
to set a definite limit on the maximum rate at which events can be processed. 

9.1. Moving the weighting and summing operations into dedicated logic 

[0125] This limitation can be avoided by placing the multiplication (weighting) and 
addition (summing) steps into dedicated logic (typically the same gate array that is 

10 performing the digital filtering and delaying operations), where they can be carried out in 
parallel once per clock cycle. FIG. 9A shows a schematic of the resultant dedicated logic 
if the four filters 205-208 of Eqn. 23c are followed by appropriate delays 210-213 to 
implement a set of cMOS secondary signal paths that include weights as well as delays. 
The weighting coefficients are applied by the multipliers 215-218 and the summation is 

15 carried out by an adder 220. The resultant, continuously processing circuit 98 thus 

becomes a drop-in replacement for the same functional block 98 in the contingent capture 
spectrometer shown earlier in FIG. 6A. A circuit of this type can be readily implemented 
in many of the higher performance field progranmiable gate arrays that are commercially 
available, for example in the Xilinx Vertex series. 

20 9.2. Order of operations is not important 

[0126] In the schematic in FIG. 9A, each secondary signal path connected to the adder 
provides an individual filter term in Eqn, 23a, appropriately delayed and weighted. Since 
the operations are independent for each path, it is clear that the order of operations may 
be rearranged in any given path without affecting the final result Ogx- Thus it is not 
25 important whether the signal is filtered, then delayed and then weighted, as the operations 
were conceptually presented, or whether some other order is adopted (e.g. weighting, 
delaying and then filtering). In practice, therefore, the operations should be arranged in 
the order that is easiest to implement. 

[0127] In the schematic in FIG. 9A, each secondary signal path connected to the adder 
30 provides an individual filter term in Eqn. 23a, appropriately delayed and weighted. Since 
the operations are independent for each path, it is clear that the order of operations may 
be rearranged in any given path without affecting the final result agx- Thus it is not 
important whether the signal is filtered, then delayed and then weighted, as the operations 
were conceptually presented, or whether some other order is adopted (e.g. weighting, 
35 delaying and then filtering). In practice, therefore, the operations should be arranged in 
the order that is easiest to implement. 
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[0128] The schematics shown in FIGS. 9B through 9E elaborate upon this point. FIG. 
9B shows a simpHfied version of the processing block 98 of FIG. 9A, where the ADC 
input signal is directed into only three secondary signal paths. Within possible processing 
operations that can be applied to the secondary signals are indicated by the sub-blocks 
5 labeled a for filtering, D for delaying, M for multiplying or weighting, S for summing, 
and C for capturing. In FIG. 9B the order of operations is the same as in FIG. 9A and 
the same in all signal paths, namely filtering, delaying, weighting, summing and 
capturing. This is a continuously processing filter, so that the capturing step requires a 
control, as we shall discuss, in order to properly obtain the energy of the input event. We 
10 note that the delays in the three paths introduce time delays of ti, t2, and ts respectively. 
In FIG. 9C the sequences of filtering, delaying and weighting operations have been 
reordered in the three secondary signal paths. As discussed above, the continuous signal 
produced by the summing operating is unchanged from FIG. 9B. 

[0129] In the circuit of FIG. 9D an additional summation point has been introduced in 
15 order to allow shorter (and typically cheaper) delays to be used. The assumption in 

moving from FIG. 9B to FIG. 9D is that the secondary signal paths were ordered so that 
the delay times ti initially satisfied the inequality ti > t2 > ts. FIG. 9D demonstrates that 
multiple summation points are allowed and that they can be distributed within the 
secondary signal paths. The circuit in FIG. 9D is also a continuously processing circuit. 

20 [0130] FIG. 9E shows an asynchronous capture circuit modeled on the circuit of FIG. 
6. Here the delays of FIG. 9B have been replaced by timed capture operations, where the 
times between captures are adjusted so that the captured secondary signal samples bear 
the same time relationships to one another as the continuous secondary signals entering 
the summation point in FIG. 9B. We observe that the summation point in FIG. 9E 

25 performs asynchronous summations because the terms that it adds arrive at different time. 
It is accordingly labeled 2a to emphasize this point. The final, single value produced by 
all four circuits of FIG. 9B - FIG. 9E after all summing and capturing operations is 
identical and this point serves to emphasize the flexibility of the invention method. 

9.3. Output of optimized TCWFS-filtering is a flat-topped trapezoid 

30 [0131] An obvious question is: what does the output of an optimized, continuously- 
filtering block look like at the summation point 220? The answer lies in our 
demonstration that the recovered values of Qgj are independent of the placement or shape 
of the step-like pulse in our filter gap, provided that the pulse's risetime is completely 
contained within the gap. This implies that the output of the filter block will be a constant 

35 over the period that the pulse's leading edge passes through the risetime sampling filter 
Go. On the other hand, since there is no output when the leading edge is not present in 



38 



any filter, the overall shape must be essentially trapezoidal, with a rising side, a flat-top, 
and a falling side. 

[0132] FIG. 10 shows modeled results demonstrating this conclusion. The figure 
shows an input pulse 232 with a 500 ns risetime. Although it is not visible on this 
5 timescale, the pulse data, and all filter data, have been decimated by a factor of 2, so the 
time step between samples is 50 ns. The summed filter output 233 is the output of a filter 
comprising three terms: a leading running average filter of length 4 ^is, a gap running 
average filter of length 1 ^is, and a lagging running average filter of length 4 ^xs. The 
baseline of this filter is exactly zero and at its peak it has 5 values (over 250 ns) that are 

10 exactly 1.000000 to one part in a million. Its rise and fall times are 4 fis. The weighted 
outputs of the leading 235, lagging 236 and gap 237 filters are also shown. These, then, 
are the cMOS secondary signals as they arrive at the summation point 220. The lagging 
filter output is shown inverted, so that it may be directly compared to the leading filter 
output. They deviate substantially only in the region where the summed filter output has 

15 its trapezoidal peak. We observe that the weighted gap filter output 237 is a small 

correction signal whose amplitude is only about 1% of the other two terms. It has been 
amplified by a factor of 20 in the figure for visibility. 

[0133] It is important to note that, while the output filter capture region 238 is easy to 
detect from the peak of the summed filter output, there are no characteristic features in 

20 any of the secondary signals (leading, lagging and gap filter outputs) upon which a 

capture might be reliably based. It is only the sum of these three signals that produces the 
desired accurate flat-top output signal. However, because filtering, delaying, and 
weighting are all deterministic processes, the time to the filter capture region following 
the arrival of a step in the input signal is well determined. Hence timing methods work 

25 very well when it is desirable to capture the individual secondary signals in contingent 
weighted sum filtering, and we have already described several circuits employing this 
principle. 

[0134] Next we will examine the behavior of a TCWFS-filter's flat-top region in more 
detail. In FIG. IIA we show the leading edges of three modeled (i.e. noise-free) input 
30 signal pulses with risetimes of 25, 250 and 500 ns (traces 240, 241, and 242 respectively), 
and a single 50 ^s decay time. The pulses were then filtered using a TCWFS-filter 
similar to the filter used to produce FIG. 10, but having a 2.80 ^is risetime and 700 ns 
gap. 

[0135] FIG. IIB shows the filter's outputs, which are trapezoidal for all three signal 
35 risetimes. FIG. IIC shows an enlargement of the trapezoids' flat-tops. Outputs 245, 

246, and 247 correspond to the inputs 240, 241, and 242 respectively). Three features are 
clear. First, the longer the risetime, the more rounded the output pulse is towards its peak. 
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Second, the length S of the flat-top is, as we would expect, equal to the difference 
between the gap length and the pulse risetime. Finally, the outputs all converge to exactly 
the same value. These flat-tops are, in fact, very flat indeed, as we can find by reference 
to Table 1. As we saw there, a high quality TCWFS-filter output may be insensitive to 
5 risetime variations at levels as low a 1 part per million, which is also a measure of how 
flat the top will be as either the risetime changes or the pulse moves about within the gap. 
Given this flatness, we can therefore capture a sample from which to recover Qgx in 
either of two ways: first by capturing a value in response to detecting a step-like pulse; 
or, second, by looking for peaks in the sunwned signal output and capturing a peak 

10 maximum. In the simplest of cases one could even do as was done in the early days of 
analog processing and simply look for output pulses with a peak detect and capture 
circuit. Clearly, in cases where there is a significant data rate, it will generally be 
preferable to also implement fast channel pulse detection and pileup inspection circuitry, 
but doing so is not a necessary part of the technique. Also, as noted earlier, capturing a 

15 "peak sample" in response to detecting a pulse may be carried out using either continuous 
or contingent TCWFS-filtering. In the latter case, as we have described, the sampling 
will actually be implemented in the cMOS secondary signal paths and these samples will 
then be combined via a "contingent" summation operation to produce the desired single 
peak value. We thereby understand that, in the contingent approach, the order of the 

20 sunmiing and capturing operations are reversed, compared to the continuous approach. In 
certain cases, the two approaches could even be mixed, with capture occurring before the 
summation operation in some signal paths and after it for others. All that is really 
required is that each secondary signal be captured after the filtering and delaying 
operations in its signal path, realizing, as noted above, that the capture method can be part 

25 of the method used to implement a required delay. 

9.4. Example implementations 
[0136] FIG. 12 shows the circuit of FIG. 6A wherein the block 98 has been replaced 
with a continuously operating TCWFS-filtering unit similar to the one shown in FIG. 9A. 
We notice that the rest of the circuit is essentially identical, with one exception. Because 

30 there is only a single value to capture, all of the digital registers (105, 109, 116, 118, and 
119) have been replaced by a single gated peak capture register 225. This register is 
triggered by the pileup test and timer module 103 in a manner exactly similar to the 
described operation of FIG. 6A in order to capture a value from TCWFS-filtering unit 98 
in the flat-top of its output. 

35 [0137] This implementation illustrates some of the flexibility that is possible in creating 
the set of secondary signal paths. In the first case, a single digital filter (021) has two 
different delay and multiply chains applied to its output in order to produce two of the 
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required set of four secondary-signals. In the second case, the gap filter (Qg) is preceded 
by its delay and weighting stages because fewer resources are required to implement them 
for the 14 bit ADC signal than for the 20 bit output from the gap filter. 

[0138] FIG. 13 shows an even simpler circuit that implements the same filter that was 
5 used to produce FIGS. 10 and 11. Its TCWFS-filtering unit 98 has only two digital 

filters: a running sum aL of length L (the peaking time) and a running sum oq of length 
G (the gap). Gl is used in two signal paths, with one path delayed by L+G with respect to 
the other to form a trapezoidal filter of risetime L (the weighting factor Ml is negative). 
gg is used in only a single path, with a delay of L, to capture the region in the other 
10 filter's gap. With appropriate weighting coefficients Ml and Mq, the summation will be 
proportional to Qgj. Note that, by scaling the weighting coefficients so that one is equal 
to unity, we avoid having to implement a third multiplication, saving resources, at the 
cost of having to rescale our spectrum later in the DSP 95 after data collection is 
complete. 

15 [0139] To demonstrate how it is done, FIG. 13 implements a different method of 
sampling the time correlated, weighted filter sum-filtering unit 98 than FIG. 12 does. 
Here the gated peak capture circuit 225 has been replaced by a gated peak detect and 
capture circuit 260 whose function is to compare each input value with a test value and, if 
the test value is smaller, to replace it with the input test, so long as the gate is asserted. 

20 Reasserting the gate sets the initial test value to zero. In this circuit the pileup test and 
timer 103 determines whether a filtered pulse is piled up or not and, if not, gates on the 
peak detect and capture circuit 260 to capture a value from the flat-top in the TCWFS- 
filtering unit's 98 output. As noted earlier, the pileup inspection part of this circuit is not 
a necessity, but has become standard in modem spectrometers because of the improved 

25 spectral quality that arises from eliminating pileup events. 

9.5, Using implicit filter delays to replace explicit delays 
[0140] As noted earlier, filtering is a convolution operation that encodes signal data 
over a time period that may either be of finite, as in finite impulse response (FIR) filters, 
or infinite duration (as in infinite impulse response (IIR) filters. By utilizing this 
30 characteristic, it is thereby possible to implement TCWFS-filters without the use of 

exphcit delay elements. The following example shows how this would be done for the 
simple running average filter used to produce FIGS. 10 and 11. That filter can be written 
as: 

+ g - (47) 

35 where 2+, Sg and 2. are the leading, gap, and lagging running average sums, respectively. 
Explicitly, these sums are: 
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i i-L i—L—M 

where the outputs of the Sg and S. filters must be delayed by L and L+M, respectively, in 
order to be combined with the output of E+. If we now define two new filters, Z+g and 
by: 

j=i~L-M+l ^ ^ j=/-2L~M+l ^ 

then it is clear that we can re-express Eqn, 47 as: 

Fj = (A-B-C) 2+ + (B-C) Z+g + C Z+gM, (50) 

where the outputs of the new filters Z+g and Z+gM can be captured simultaneously with 
the output of filter Z+. Thus, by using the implicit delays in the filters Z+g and Z+gM, we 
10 have eliminated the need for the explicit time delays required by the formulation of Eqn. 
47. 

10. Two-Pole example application to CsKTO scintillators 

[0141] In the majority of applications described above, the signals were nominally 
single-pole outputs with perhaps a small second-pole term as an error term in the desired 

15 output. The TCWFS-filtering method was originally developed to deal with such signals, 
and it does so well. However, the method is not restricted to such signals, as we noted 
earlier and the following example demonstrates. In this case we will look at signals from 
CsI(Tl) scintillators coupled either to photodiodes or photomultiplier tubes. It has been 
known for some time that CsI(Tl) is one of a class of scintillators whose light output 

20 decays with at least two different decay time components [SKULSKI - 2001], a short 
component of about 600 ns and a long component of about 4 ^.s. FIG. 14 shows an 
output pulse from CsI(Tl) when stimulated by a gamma-ray. The shorter and longer 
components are clearly visible in the output, which is readily described as a two-pole (2- 
P) device. As a scintillator, CsI(Tl) has several properties that make it superior to the 

25 more commonly employed Nal(Tl) scintillator, which has only a single decay time. The 
presence of the second decay time, however, has so far precluded the useful application of 
CsI(Tl) at significant counting rates because it interferes sufficiently with the operation of 
existing digital and analog shaping filters, so that energy resolution degrades rapidly 
above 10,000 counts/sec (cps) and becomes unusable above 20,000 cps. 

30 [0142] To work with CsI(Tl) at higher rates we developed a TCWFS-filter comprising 
five running average terms. The light output from CsI(Tl), as seen in FIG. 14, is clearly a 
step-like pulse. However, as discussed above, here it is the area under the curve (the total 
number of photons emitted) that represents the energy deposited in the scintillator, rather 
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than the pulse amplitude as in other common detectors. As a result, the derivation of the 
weighting functions in the TCWFS-filter must be modified slightly from the derivation 
presented in §4. We therefore present it in compact form below. 

[0143] We model the light output as a two-pole device, with two decay times and T2 
5 and two decay coefficients b^ and b2, per §4. The amplitudes Q^^ and Q23 represent the 
residual amplitudes of these components from any preceding absorption events at time 3, 
following the gap period g. Qgi3 and Qg23 represent the increases in the amplitudes of 
these terms at the same time instant, is the light emitted during the signal's risetime in 
the gap period g and equals the sum of the light emitted by the two components, Og^ and 
10 ag2. For simplicity, only three filter lengths, L, M, and g, are used to implement the five 
running sum filters, as shown in FIG, 15. Given that the measured light output at time i 
is y., the outputs of the five running sum filters Oj are then: 

= bf(L-.M^g) Ail + Q23 b2-<L+M-^S) (51a) 

<^2 = Ql3 bf^^^^ \m + Q23 V^^^^^ (51b) 
15 03 = bf(g) + Q23 b2-(8> A2g + CTg (51c) 

04 = (Qi3 + Qgia) ^iL + (Q23 + Qg23)A2L (5 Id) 

<^5 = (Ql3 + Qgl3) Ail + (Q23 + Qg23) t>2^ (51e) 
where the Ay terms have the same definitions as in § 4. 

[0144] If we define = a^, q2 = Qg^, q3 = Qg23' ^4 = Ql3' ^5 = ^23' 
20 similarly to Eqns. 23a and 23b, we can write Eqn. 51 in matrix form and solve as: 

a = J-q, (52a) 
so that q = J"* • a, (52b) 

[0145] Now, the total emitted light, taken as a measure of the deposited energy Qg^, can 
be expressed in the model parameters Q, where r^ and r2 have the same definitions as in 
25 § 4, as: 

QgT = + ^1 Qgl3 + ^2 Qg23- <53) 
What Eqn. 53 represents is our best estimate, using only filters of length g, M, and L, of 
the light that we would have obtained if we had actually integrated to infinity, rather than 
only over the finite time represented the filter lengths. We can replace the Q*s from Eqn. 
30 52 and collect terms to obtain: 

V = ^'U^i-"-! 2-^2"/^i ^^2 h^i'^r or (54a) 
* i = l 1=1 1 = 1 
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i-\ 1=1 
where the values are the weighting coefficients to be applied to the 5 time correlated, 
filtered secondary signals that we derive from our three primary filters by appropriate 
delays. Thus is a running average of length L, delayed by L+2M+g; O2 is a running 
5 average of length M, delayed by L+M+g; 03 is a running average of length g, delayed by 
L+M; O4 is a running average of length M, delayed by L; and 05 is a running average of 
length L with no delay. 

[0146] The lengths L = 1.6 ^s, M = 0.8 |lls, and g = 0.6 jis were found to work well. For 
example, 1.6 [is filter length (the "peaking time") was varied to see what its effect was on 

10 energy resolution and high rate performance. In all cases the correct weighting functions 
were computed for the filter lengths and CsI(Tl) decay times using Eqns. 51 - 54. Results 
of energy resolution versus peaking time are shown in FIG. 16 for a 1" cubic scintillator 
exposed to ^^^Cs radiation at 662 keV. The energy resolutions obtained at 19 kcps are 
essentially identical to values obtained for a crystal of this size with traditional processing 

15 electronics. In contrast to earlier filtering methods, however, the energy resolution now 
degrades quite slowly with rate. At 3.2 peaking time, resolution increases only 10% 
(6.0 to 6.6%) going from 19 kcps to 104 kcps and only 7% (9.1 to 9.7%) at 1.2 ^is 
peaking time. Slightly larger increases are found going to 207 kcps. These energy 
resolutions are still useable, however, and essentially extend the count rate capability of 

20 CsI(Tl) by over an order of magnitude compared to earlier filtering methods. 

[0147] Finally, it is also known that the amount of light CsI(Tl) emits in its two decay 
components varies depending upon the identity of the particle depositing energy in the 
crystal. FIG. 17 shows typical traces from a gamma-ray, a proton, and an alpha particle 
depositing the same amount of energy in CsI(Tl). The change in curve shape with 

25 particle identity allows "Particle ID" or FID schemes to be implemented, as described in 
[SKULSKI -2001] and elsewhere. A typical method is to compare the area under the 
first part of the curve (e.g. the first 1.5 \Jis) to the total area under the curve (e.g. its 
energy). Obviously, if the energy cannot be determined accurately at higher counting 
rates, because of the presence of two decay terms as discussed above, then the FID will 

30 not be accurate at higher counting rates either. The situation can be improved, therefore, 
by using an energy value determined by the TCWFS-filtering methods described herein. 
[0148] In the process of solving for Qg^, the TCWFS-filtering method also produces 
terms that can be directly used to implement FID. As a first example, if the ratio of light 
in the two decay modes is characteristic of the stimulating particle, then Qgi3/Qg23 

35 direcdy reflects this ratio and can be used as a FID parameter. Altematively, since 
inverting matrix J introduces some correlation between Qgi3 and Qg23» can also 
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compare the area under the early part of the light output curve to the total area. In this 
case CTg/Qgj can be used as the PID parameter. The fourth filter, a^, also measures early 
light output, but it rides on the light output of any earlier events when operating at high 
count rates. However, once the vector Q has been found, it can then be corrected for this 
error by subtracting these terms to get the true area, A4, under arising from the new 
event, from Eqn. 51d: 

A4 = 04 - Ail - Q23 (55) 
allowing A4/Qgj to be used as the PID parameter. Thus the TCWFS-filtering method 
also allows PID techniques to be extended to high counting rate situations. 
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12. Conclusion 

[0150] The foregoing description of specific embodiments has been presented for 
purposes of illustration and description. It is not intended to be exhaustive or to limit the 
5 invention to the precise forms described, and obviously, many modifications and 

variations are possible in light of the above teaching. These embodiments were chosen 
and described in order to best explain the principles of the invention and its practical 
application to thereby enable others in the art to best utilize the invention in various 
embodiments and with such modifications as best suit the invention to the particular uses 
10 contemplated. 
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[0151] As a first example, while all the described embodiments used trapezoidal, 
triangular, or running average filters, there is no such restriction in the technique. As 
shown, any filter describable by a set of weighting constants applied to successive signal 
samples, or its analog equivalent, can be used. The method will therefore work with 
5 essentially any realizable filter. Second, while the specific embodiments employed two, 
three, or four filters, there is no specific limit on the number of filters which can be used. 
Even a single filter can be used if its output is correctly sampled at several times. If the 
number of filters N exceeds the number of base amplitudes M, then matrix inversion 
methods can be used to find least squares best estimates of the M values, rather than exact 

10 values. Third, while all the filters used in the specific embodiments were symmetric, this 
is not a requirement of the method, as may be seen from the derivations. Fourth, while 
specific configurations of electronic circuits were shown in the implementations, other 
configurations can readily be devised which accomplish the same functions. The number 
of ADC bits or value lengths captured to the registers could be varied, for example. Fifth, 

15 even the use of electronic circuits is not required to implement the method. As per the 
initial derivations, for example, if successive ADC values are instead read into computer 
memory, all the subsequent processing steps can be carried out in software. Sixth, while 
all of the described embodiments used one or more explicit delay elements to create the 
desired time correlations between filter outputs, as we showed in § 9.5, it is not difficult 

20 to produce designs where delays implicit in the filters themselves can be made to serve 
this function. Seventh, in some cases an unfiltered sample of the primary signal, i.e. a 
unity filtered sample, may be one of the secondary-signals. Eighth, the method can also 
be applied to signals which are formally equivalent to outputs from multipole filters, 
whether the signal source is actually a single multi-pole filter or not. Superconducting 

25 bolometers and certain classes of scintillators produce such output signals, for example. 
[0152] Therefore, the above description should not be taken as limiting the scope of the 
invention, as defined by the appended claims. 
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